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Introductory  Note 

Research  has  been  a  continuing  part  of  the  programs  of 
the  Naval  Weapons  Center  from  its  earliest  days.  Today  the 
Center's  mission,  involving  air  weapons  and  tactical  missiles, 
requires  the  Center  to  maintain  a  broad  program  of  research  in 
physics,  chemistry,  mathematics,  computer  science,  and  the 
engineering  sciences  of  detonation  physics  and  aerother- 
mochemistry.  The  articles  in  this  issue  of  the  Naval  Research 
Reviews  provide  examples  and  summaries  of  some  of  the 
Center’s  research  work. 

Research  efforts  of  the  Naval  Weapons  Center  are  inte¬ 
grated  with  the  Center's  systems  development  activities. 
Thus,  it  is  useful  to  view  the  program  in  terms  of  potential 
applications  of  its  results. 

The  research  program  in  physics  is  focused  on  optical, 
microwave,  and  electronic  devices  and  materials.  Applica¬ 
tions  of  these  efforts  are  in  targeting,  missile  guidance,  fuzing, 
and  laser  systems. 

Research  in  chemistry  encompasses  a  broad  spectrum  of 
applications,  including  optical  and  electronic  materials,  explo¬ 
sives.  propellants,  polymers,  electrochemical  devices  for 
energy  conversion  applications,  and  instrumental  analysis  in 
support  of  all  of  these  efforts 

Applications  of  efforts  in  mathematics  and  computer  sci¬ 
ence  include  systems  modeling  and  analysis,  signal  proces¬ 
sing.  control  theory,  stochastic  processes,  and  applied  artifi¬ 
cial  intelligence 

Explosives,  warheads,  and  terminal  effects  are  the  sub¬ 
jects  of  interest  in  the  research  program  in  detonation  physics, 
while  missile  propellants  and  propulsion  systems  are  of  in¬ 
terest  in  the  aerothermochemistry  research  program. 

Research  scientists  and  engineers  arc  generally  given  the 
freedom  to  explore  research  topics  of  their  own  choosing 
within  this  matrix  of  subjects,  applications,  and  disciplines. 
While  the  Center's  research  program  is  not  large,  it  produces 
annually  some  SO  to  60  refereed  publications  and  a  similar 
number  of  presentations.  Funding  support  for  the  program 
comes  primarily  from  the  Naval  Material  Command,  the  Of¬ 
fice  of  Naval  Research,  the  naval  systems  commands,  and  the 
Defense  Advanced  Research  Projects  Agency. 

Dr.  Edwin  B.  Royce 
Head,  Research  Department 
Naval  Weapons  Center 


2 


Naval  frtwn*  Rrvie hi 


Tectonics, 
Seismicity,  and 
Yolcanism 
at  The  Naval 
Weapons 
Center 

by 

Glenn  R.  Roquemore 
and 

John  T.  Zellmer 
Naval  Weapons  Center 


Introduction 

The  Naval  Weapons  Center  (NWC)  at  China  Lake  oc¬ 
cupies  a  part  of  California  known  for  its  young  geologic  fea¬ 
tures  and  high  level  of  seismic  activity.  The  Coso  Range  and 
the  Indian  Wells  Valley,  which  lie  within  the  boundaries  of 
NWC.  are  located  in  the  southwestern  comer  of  the  Basin  and 
Range  physiographic  province.  The  Coso  Range  encompasses 
1 ,920  square  kilometers;  to  the  west  is  the  Sierra  Nevada,  and 
to  the  east  is  the  Argus  Range.  Owens  Lake  terminates  the 
Coso  Range  to  the  north.  The  Indian  Wells  Valley  covers 
1 ,020  square  kilometers  to  the  south  of  the  Coso  Range  (Figure 
I).  The  largest  population  center  and  most  of  NWC’s  labora¬ 
tory  facilities  are  located  in  the  Indian  Wells  Valley.  This 
article  describes  the  tectonics,  seismicity,  and  volcanism  of  the 
Coso  Range  and  the  Indian  Wells  Valley,  and  the  recent  im¬ 
pacts  that  faulting,  crustal  deformation .  and  volcanic  activity 
have  had  on  NWC. 
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The  Coso  Range 
has  a  history  of 
volcanism  and 
earthquakes. 


Figure  1 

Index  map  of  Coso  Range  and  Indian  Wells  Valley,  showing 
NWC,  locations  of  volcanos,  and  epicenters  of  1 946  and  1 96 1 
earthquakes. 
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The  Coso  Range 

The  Coso  Range  has  a  history  of  volcanism  and  earth¬ 
quakes.  The  eruption  of  high-silica  rhyolite  occurred  at  least 
38  times  over  the  last  million  years,  and  basaltic  eruptions 
ceased  only  30,000  years  ago.  Duffield,  et  al.\  have  shown 
that,  in  all,  about  33  cubic  kilometers  of  lava  have  erupted 
during  the  last  4  million  years.  The  Coso  Range  has  been 
broken  by  Basin  and  Range  extensional  faulting,  as  seen  in  the 
extensive  horst  and  graben  development  and  right-slip  dis¬ 
placement.  Bacon2  has  shown  that  the  eruption  of  both  basalt 
and  rhyolite  in  the  Coso  Range  has  been  controlled  by  regional 
tectonic  extension. 

Roquemore3  has  shown  that  earthquakes  with  mag¬ 
nitudes  in  the  M  =  8.0  range  may  have  occurred  in  the  Coso 
Range  within  the  last  1,600  years,  and  0.3  meters  of  normal 
displacement  along  the  Airport  Lake  fault  probably  occurred 
in  1872.  Walter  and  Weaver4  have  reported  that  4,800  earth¬ 
quakes  of  M  =  1.0  to  M  =  3.9  occurred  within  the  Coso 
volcanic  field  from  1975  to  1977.  Earthquakes  of  M  a  5.0 
have  occurred  peripheral  to  the  Coso  Range.  Some  of  these 
earthquakes  include  the  1872  Owens  Valley  earthquake  (M  = 
8.3)  and  the  1946  Walker’s  Pass  earthquake  (M  =  6.3).  The 
majority  of  local  earthquakes  occurring  within  the  Indian 
Wells  Valley  have  been  Ms  5.2. 

According  to  Combs5 ,  the  Coso  Range  contains  abundant 
fumaroles,  and  high  heat  flow  is  associated  with  them.  Tele- 
seismic  P-wave  delays  centered  near  the  heat  flow  anomaly 
indicate  a  potential  for  partially  melted  material  at  a  depth  of 
between  8  and  20  kilometers.  According  to  Reasenberg,  et 
al.b,  an  intense,  low-velocity  body  is  located  under  Devils 
Kitchen  (Figure  2).  This  low-velocity  body  is  approximately  5 
kilometers  wide  on  top  and  becomes  increasingly  elongated  in 
the  north-south,  or  N25°  E  direction  with  depth.  A  maximum 
velocity  contrast  between  the  partially  melted  material  and  the 
surrounding  rock  is  attained  between  10  and  17.5  kilometers  in 
depth.  The  heat  source  in  the  Coso  volcanic  field  appears  to  be 
centered  near  Devils  Kitchen. 

Most  of  the  dozens  of  known  fumaroles  in  the  Coso 
Range  are  associated  with  faults.  Faults  in  the  Coso  Range  that 
are  associated  with  hot  springs  have  significant  components  of 
normal  slip.  In  general,  these  faults  are  segments  of  the  Airport 
Lake  fault  zone  and  contain  the  highest  concentration  of  hot 
springs  in  the  region.  Both  the  Airport  Lake  and  the  right-slip 
Little  Lake  faults  have  displacements  and  strikes  consistent 
with  the  regional  tectonic  stress  pattern. 


Hgun  2 

Location  of  low-velocity  body  under  Devils  Kitchen. 


Figure  3 


Aerial  view  ot  the  bttle  Lake  fault  zone,  showing  typical 
geomorphic  features  and  offset  in  a  400.000-year  old  basalt 
flow. 
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Figure  4 

Locations  of  Airport  Lake  fault  zone,  Uttle  Lake  fault  zone,  and 
White  Hills  anticline. 
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Little  Lake  and  Airport  Lake  Faults 

The  Little  Lake  and  Airport  Lake  faults  are  the  two  major 
active  faults  in  the  Indian  Wells  Valley  and  the  Coso  Range. 
These  faults  are  components  of  the  regional  right-slip  shear 
and  the  east-west  extension  that  characterize  the  tectonics  of 
the  western  Basin  and  Range  physiographic  province.  The 
Little  Lake  and  Airport  Lake  fault  zones  developed  in  response 
to  the  same  regional  tectonic  stress  field.  Although  these  fault 
zones  share  a  common  origin,  they  represent  different  styles  of 
deformation.  The  Little  Lake  fault  is  predominantly  right  slip, 
but  develops  a  normal-slip  component  along  its  central  and 
southern  segments.  It  also  has  many  geomorphic  features  simi¬ 
lar  to  those  of  the  San  Andreas  fault  located  on  the  western 
margin  of  California.  The  Airport  Lake  fault  displays  normal 
slip  displacement  and  contains  numerous  nested  grabens. 

From  the  northwest  comer  of  the  Indian  Wells  Valley,  the 
Little  Lake  fault  splays  southeast  from  the  main  trace  of  the 
Sierra  Nevada  frontal  fault  and  trends  diagonally  across  the 
valley  to  at  least  the  NWC-Ridgecrest  border.  A  seismic  re¬ 
fraction  study  by  Zbur7  suggests  that  the  bedrock  expression  of 
the  Little  Lake  fault  continues  southward  to  the  Garlock  fault. 
The  fault  pattern  is  well-exposed  on  a  basalt  flow  in  the  north¬ 
west  part  of  the  Indian  Wells  Valley  (Figure  3).  Rhombic 
depressions,  pressure  ridges,  and  linear  troughs  are  typical 
along  this  portion  of  the  fault.  As  the  fault  enters  the  deep 
alluvial  basin  of  the  Indian  Wells  Valley,  its  expression 
changes  to  short,  segmented  scarp  traces  that  form  an  en  eche¬ 
lon  pattern. 

The  Airport  Lake  fault  strikes  north-south  through  the 
Coso  Range  and  the  Indian  Wells  Valley  for  over  35  kilome¬ 
ters  (Figure  4).  The  northern  sigment  of  the  fault  is  a  narrow 
zone  less  than  0.5  kilometers  wide.  The  central  segment 
widens  to  8  kilometers,  forms  Coso  Basin,  and  extends  into  the 
Indian  Wells  Valley.  Cinder  cones,  basaltic  dikes,  and  sym¬ 
metrically  stepped  grabens  occur  within  the  fault  zone  (Figure 
5).  The  fault  pattern  bounding  the  zone  is  in  places  suggestive 
of  oblique  slip;  however,  within  the  fault  zone  the  patterns  are 
clearly  normal  slip  and  are  associated  with  graben  formation. 
Vertical  displacement  near  Coso  Valley  totals  at  least  600 
meters  within  a  zone  5  kilometers  wide. 

The  Little  Lake  and  the  Airport  Lake  faults  intersect  near 
the  center  of  the  Indian  Wells  Valley  (Figure  4).  This  intersec¬ 
tion  is  the  locus  of  repeated  seismic  swarms  that  have  been 
occurring  since  April  1981.  Because  of  the  high  level  and 
characteristics  of  seismicity  and  the  notable  changes  in  the 
surface  expression  of  the  faults,  the  intersection  area  has  be¬ 
come  the  subject  of  much  study.  These  changes  suggest  that 
much  of  the  regional  extension  associated  with  the  Airport 
Lake  fault  is  being  transferred  to  the  Little  Lake  fault  as  right- 
slip  displacement,  both  within  and  to  the  south  of  the  zone  of 
intersection. 
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Geodetic  and  Level-Line  Surveys 

Survey  data  obtained  at  NWC’s  supersonic  naval 
ordnance  research  track  and  within  the  geodetic  and  level-line 
trilateration  network  provide  information  about  the  regional 
and  local  tectonic  stress  fields  and  the  related  deformation 
occurring  in  the  Coso  Range  and  the  Indian  Wells  Valley  area. 
These  data,  in  combination  with  tectonic  studies  such  as  those 
done  by  Wright*,  indicate  a  regional  stress  field  (Figure  4) 
resulting  in  right-slip  shearing  (evidenced  in  the  Little  Lake 
fault),  extensional  deformation  (evidenced  in  the  Airport  Lake 
fault),  and  folding  (evidenced  in  the  White  Hills  anticline). 
Similar  stress-related  features  on  a  smaller  scale  are  also  found 
in  the  Coso  Range,  in  the  Indian  Wells  Valley,  and  along  the 
research  track. 

Periodic  and  highly  precise  surveys  of  the  research  track 
are  used  to  maintain  track  alignment  tolerances  for  supersonic 
rocket  sleds.  Surveys  of  the  7-kilometer-long  research  track 
and  of  Coso  Range  show  that  the  region  is  actively  deforming. 
Surveys  suggest  that  between  1977  and  1978,  the  track  up¬ 
lifted  29  millimeters  at  its  northern  end  and  also  may  have 


Figure  3 

Aerial  view  of  the  Airport  Lake  fault  zone,  showing  a  typical 
Graben. 


experienced  minor  lateral  deformation  and  shortening.  The 
track  is  being  resurveyed  to  determine  if  additional  deforma¬ 
tions  have  occurred. 

Repeated  surveys  in  the  Coso  Range  have  included  both 
leveling  to  detect  elevation  changes  and  trilateration  to  detect 
horizontal-slip  displacements.  On  13  March  1982  a  possible 
spasmodic  tremor  (indicating  volcanically  induced  seismic  ac¬ 
tivity)  was  recognized  near  Cactus  Peak  in  the  Coso  Range.  To 
date,  it  is  not  known  if  the  tremor  actually  occurred;  however, 
several  existing  geodetic  and  level  lines  were  resurveyed  and 
new  survey  lines  were  added.  Leveling  surveys  performed 
along  a  20-kilometer-long  line  across  the  Coso  volcanic  field 
indicate  that  a  30-millimeter  uplift  occurred  between  1975  and 
1982.  The  result  of  the  trilateration  network  surveys  indicate 
that  a  slight  north-south  shortening  and  east-west  lengthening 
may  have  occurred  between  1975  and  1982.  One  hypothesis  to 
explain  the  uplift  and  geodetic  deformation  along  the  track  is 
that  a  sill  of  magma  may  have  intruded  into  the  upper  crust. 
This  hypothesis,  proposed  by  Savage’,  is  being  tested  by  con¬ 
tinuing  measurements  and  detailed  studies  of  the  earthquake 
data. 

Because  of  the  scientific  interest  generated  by  the  earth¬ 
quake  swarms  occurring  since  April  1981,  the  trilateration 
network  was  recently  expanded  to  include  45  lines  having  a 
combined  length  of  about  441  kilometers.  The  total  area  cov¬ 
ered  by  the  network  is  about  230  kilometers. 

Seismicity 

The  earthquake  sequences  that  began  during  April  1981 
marked  the  end  of  a  20-year  period  of  seismic  quiescence  at 
NWC.  The  Indian  Wells  Valley  area  had  not  experienced  any 
earthquakes  of  M  2  4  since  1961.  During  the  period  from 
April  1981  to  October  1982,  there  were  five  earthquake  se¬ 
quences  with  maximum  magnitudes  ranging  from  M  =  4.0  to 
M  =  5.2  and  several  other  sequences  with  earthquakes  of 
lesser  magnitude .  These  sequences  contained  1 ,029  detectable 
earthquakes,  169  of  which  exceeded  M  =  2.5  and  1 1  of  which 
equaled  or  exceeded  M  =  4.0.  The  seismic  record  for  the 
period  of  1  December  1980  to  5  October  1982  is  shown  in 
Figure  6.  Although  the  data  set  is  complete  only  through  5 
October  1982,  a  high  number  of  low-magnitude  earthquakes 
are  still  occurring  as  of  December  1982. 

The  seismic  record  and  the  epicentral  regions  shown  in 
Figure  4  indicate  three  general  trends  associated  with  the 
earthquake  sequences.  With  each  successive  sequence,  the 
tendency  has  been  for  the  total  number  of  earthquakes  to  in¬ 
crease,  for  the  maximum  magnitude  to  increase,  and  for  the 
epicentral  region  to  migrate  southward.  These  sequences  typi¬ 
cally  continue  for  several  days  or  weeks.  The  first  sequence, 
occurring  in  April  1 98 1 ,  contained  five  earthquakes  and  had  a 
maximum  magnitude  of  4.0.  The  second  sequence  occurred  in 
January  1982  and  contained  8  earthquakes,  the  largest  of 
which  was  M  =  4.2  The  third  sequence,  occurring  in  March 
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Figure  • 

Seismic  record  for  the  period  of  1  December  1 980  through  5 
October  1982. 


1982,  contained  25  earthquakes  and  reached  a  maximum 
magnitude  of  M  =  4.5.  The  fourth  earthquake  sequence,  oc¬ 
curring  during  the  period  of  April  to  August  1982,  did  not 
appear  to  fit  the  trend,  but  calculations  indicate  that  the  total 
amount  of  seismic  energy  released  exceeded  that  of  the  March 
1982  sequence.  Even  though  the  maximum  magnitude  did  not 
increase,  both  the  number  of  earthquakes  in  die  sequence  and 
the  energy  released  showed  an  increase.  The  fifth  and  most 
recent  sequence  began  in  late  September  1982.  This  one  con¬ 
tained  at  least  81  earthquakes  and  reached  a  maximum  mag¬ 
nitude  of  5.2.  Because  data  are  not  yet  available  for  the  period 
after  5  October,  the  total  number  of  earthquakes  associated 
with  the  sequence  and  the  duration  of  the  sequence  are  not  yet 
known. 

Along  with  the  trends  described  above,  the  earthquake 
epicenters  also  changed  locations  with  each  succeeding  se¬ 
quence.  The  epicentre!  regions  (shown  in  Figure  4)  migrated 
southward  and  showed  changes  in  orientation.  The  shapes  and 
locations  of  die  epicentral  regions  suggest  that  the  earlier  se¬ 
quences  were  associated  predominandy  with  the  Little  Lake 
fault  and  that  the  later  sequences  may  have  had  au stronger 
affinity  to  the  Airport  Lake  fault.  First-motion  studies  of  the 
earthquakes  have  so  far  been  limited,  but  suggest  that  right-slip 
motion  has  been  the  dominant  mechanism. 


Discussion 

Geodetic  and  Level- Line  Surveys 

The  styles,  magnitudes,  and  rates  of  strain  observed 
along  the  research  track  and  within  the  geodetic  and  level-line 
trilateration  net  are  indicative  of  the  deformations  occurring 
within  the  Coso  Range  and  the  Indian  Wells  Valley  region  as  a 
whole.  The  data  obtained  so  far  indicate  that  portions  of  the 
Coso  Range  and  Indian  Wells  Valley  region  are  experiencing 
north-south  compression,  east-west  extension,  and  vertical  up¬ 
lift  at  rates  significantly  higher  than  those  observed  in  much  of 
the  tectonically  active  western  United  States. 

Faulting  and  Seismicity 

Certain  areas  of  California  appear  to  be  coming  out  of  a 
20-year  lull  in  seismicity;  NWC  occupies  one  of  these  areas. 
Areas  to  the  north,  including  Mt.  Shasta,  Mt.  Lassen,  and  the 
Long  Valley  Caldera  near  Mammoth  Lakes,  all  seem  to  be 
showing  signs  of  renewed  seismic  activity.  The  Coso  range  is 
probably  the  southernmost  volcanic  field  included  in  that  activ¬ 
ity. 

The  high  level  of  seismicity  that  began  during  April  1981 
resulted  from  renewed  movement  within  the  Little  Lake  and 
Airport  Lake  fault  system.  This  movement  indicates  that  the 
build-up  of  regional  and  local  stresses  has  exceeded  the 
strength  of  the  rock  in  die  fault  zones  and  that  stress-relief 
adjustments  are  occurring.  The  majority  of  the  adjustments,  as 
shown  by  the  earthquake  epicentral  regions  in  Figure  4,  are 
located  in  the  zone  where  die  two  faults  intersect,  increased 
high  levels  of  seismic  activity  in  the  intersection  zone  will  not 
necessarily  prevent  stress  buildup  and  release  along  other 
segments  of  the  Little  Lake  and  Airport  Lake  fault  zones. 
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Higher  levels  of  seismic  activity  can  probably  be  expected  to 
occur  within  this  zone.  Roquemore  's  analysis  of  paleoseismic- 
ity  along  the  northern  segment  of  the  Little  Lake  fault  suggests 
that  earthquakes  with  magnitudes  of  M  =  8.0  or  higher  may 
have  occurred  within  the  last  1 ,600  years. 

Future  Concerns 

If  the  volcanic  field  experienced  spasmodic  tremor  and 
magmatic  intrusion  near  Cactus  Peak  in  the  Coso  Range,  crus¬ 
tal  adjustments  (in  the  forms  of  earthquakes  on  the  Little  Lake 
and  Airport  Lake  faults)  may  be  relieving  enough  of  the  stress 
in  the  fault  zone  to  allow  for  the  propagation  of  magmatic  dikes 
or  sills  toward  the  surface. 

The  recent  earthquake  sequences  occurring  in  the  Indian 
Wells  Valley  area  and  the  trends  exhibited  by  them  are  events 
that  have  not  been  observed  before  in  this  region.  The  relation¬ 
ship  of  these  earthquakes  to  future  seismic  activity  is  not 
known.  However,  these  events,  in  combination  with  the  know¬ 
ledge  of  the  regional  tectonic  stress  fields  of  the  Indian  Wells 
Valley  area,  have  caused  the  U.S.  Geological  Survey  and 
California  Office  of  Emergency  Services  to  jointly  issue  a 
warning  that  conditions  may  be  favorable  for  earthquakes  of 
larger  magnitude  (M  2:  5.2)  to  occur  within  the  year. 
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Introduction 

The  quest  for  military  weapons  of  high  accuracy  at  long 
standoff  ranges  has  led  to  extensive  research  and  development 
efforts  to  improve  ram  jet  missile  propulsion  systems  as  well  as 
other  concepts  that  utilize  the  weight-saving  fuel  efficiencies 
possible  through  the  use  of  airbreathing  combustion  systems. 
A  part  of  the  effort,  which  is  under  way  at  the  Naval  Weapons 
Center,  involves  research  into  methods  for  improving  air- 
breathing  combustion  and  thus  increasing  its  viability  for  use 
in  weapon  systems.  The  studies  described  in  this  article  are 
concerned  with  combustion  instability  in  liquid  fuel  ramjet 
combustors. 

Additional  related  research  efforts  are  under  way  that  deal 
with  turbulent  mixing  and  combustion  in  airbreathing  systems 
such  as  gas  generator  ramjets  (ducted  rockets),  solid  fuel  ram¬ 
jets.  and  projectiles  assisted  by  external  burning  drag  reduction 
or  thrust  enhancement.1  6  The  knowledge  gained  from  these 
efforts  and  from  three  decades  of  experience  in  solid  propellant 
rocket  motor  combustion  instability  form  a  background  of 
expertise  for  the  studies  described  here. 

Combustion  Instability  in  Ramjets 

It  is  generally  known  that  combustion-driven  oscillations 
can  occur  in  a  wide  variety  of  fuel-buming  devices.  Acoustic 
instability  involves  interactions  between  acoustic  wave  modes 
with  the  combustor  and  combustion  processes  in  such  a  way 
that  minor  flow  disturbances  in  the  combustor  gas  are 
amplified  into  a  standing  acoustic  wave.  Acoustic  instability  in 
propulsion  systems  can  interfere  with  missile  performance, 
making  instability  a  frequent  problem  in  the  design  of  propul¬ 
sion  systems.  Such  combustion-driven  oscillations  also  occur 
in  ramjet  engines .  In  tests  of  ramjet  engines  of  compact  design , 
for  example,  low-frequency  (100-  to  500-  Hz),  combustion- 
induced,  longitudinal  pressure  oscillations  have  been  ob¬ 
served.’*  If  these  oscillations  are  of  sufficiently  high 
amplitude,  they  can  cause  inlet  unstart  and  result  in  loss  of 
engine  performance. 

Over  the  past  30  years,  NWC  has  developed  expertise  in 
combustion  instability  in  solid  propellants.  The  goal  of  this 
research  is  to  build  a  technology  base  to  eliminate — or  at  least 


control — instability  so  that  acoustic  pressure  amplitudes  can 
be  lowered  to  acceptable  levels.  The  progress  of  the  effort  in 
the  solid  propellant  area  has  led  to  methods  for  testing  the 
acoustic  response  of  solid  propellants  and  to  the  prediction  of 
acoustic  instability  of  solid  propellant  motors.'""  Experience 
gained  from  these  studies,  along  with  applications  of  high- 
frequency  response  instrumentation  and  data  processing,  has 
been  communicated  to  members  of  the  ramjet  engine  com¬ 
munity." 

In  recent  studies  of  ramjet  combustion  instability,  ex¬ 
perimental  and  theoretical  means  were  used  to  determine  the 
relative  importance  of  major  combustor  operational  conditions 
on  pressure  oscillation  amplitudes  to  ascertain  the  acoustic- 
wave  structure,  to  determine  the  shock  oscillation  properties, 
and  to  provide  insight  into  shock  and  acoustic  wave  interac¬ 
tion.  These  experiments  were  performed  in  a  coaxial  dump 
combustor  of  simple  design.  Experiments  to  determine  pres¬ 
sure  wave  structure  in  a  more  complex  ramjet  with  side- 
mounted  inlets  (side  dump)  were  also  pursued,  but  are  not 
reported  here.1' 

Experimental  Arrangement 

The  laboratory-sized  5-in. -diameter  heavy-walled  ramjet 
combustor  used  in  these  studies  is  shown  in  Figure  1 .  Condi¬ 
tions  for  the  majority  of  testing  conducted  to  date  are  shown  in 
Figure  2. 

Experimental  Results  and  Interpretation 

Research  at  NWC  has  concentrated  on  three  major  areas: 

( 1 )  the  effect  of  ramjet  operational  conditions  on  engine  pres¬ 
sure  oscillations,14  (2)  identification  of  acoustic  modes  that  can 
be  excited,1'  and  (3)  the  response  of  the  inlet  shock  to  acoustic 
pressure  oscillations  in  the  inlet. 16  Another  area  of  interest  on 
which  work  has  been  initiated  is  interaction  of  the  acoustic 
field  with  the  coherent  vortex  structure  generated  in  the  shear 
layer  near  the  plane  of  the  dump  and  in  the  recirculation  zone, 
(see  Figure  I). 
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Effect  of  Operating  Conditions 

To  be  of  maximum  military  use,  a  ramjet  engine  must 
operate  over  a  range  of  speeds  and  altitudes.  Those  flight 
conditions  result  in  changes  of  air  mass  flow  at  the  inlet  with 
related  changes  in  air-to-fuel  ratio,  chamber  pressure,  and  dif¬ 
fuser  exit  Mach  number,  all  of  which  can  affect  combustion 
instability.  In  addition,  design  parameters  such  as  engine 
geometry,  type  of  fuel,  fuel  injection  pattern,  and  fuel  penetra¬ 
tion  into  the  inlet  air  flow  can  affect  instability.  Experiments  at 
NWC  have  been  designed  to  test  the  effects  of  those  variables. 
One  significant  finding  was  that  longitudinal  acoustic  instabil¬ 
ity  was  relatively  insensitive  to  fuel  type,  permitting  tests  to  be 
performed  with  a  gaseous  fuel  (ethylene,  C2H4)  rather  than  a 
liquid  fuel  (such  as  JP-4).  On  the  other  hand,  optimum  condi¬ 
tions  for  instability  (highest  acoustic  amplitudes)  were  rather 
narrow,  requiring  80  psia  chamber  pressure,  a  5-lb/s  inlet 
flow,  and  a  diffuser  exit  Mach  number  of  approximately  0.4. 

One  of  the  more  significant  results  in  terms  of  probable 
instability-driving  mechanisms  was  the  effect  of  fuel  injection 
velocity  (penetration  into  the  inlet  air  flow)  on  the  magnitude 
of  the  pressure  oscillations  (Figure  3).  The  injection  velocity 
affects  fuel  entrainment,  fuel  mixing  in  the  shear  layer,  and 
mixing  and  combustion  in  the  recirculation  zone  (see  Figure 
I),  which  are  of  great  importance  in  driving  the  combustion 
instability.  This  subject  is  under  extensive  investigation. 


Identification  of  Acoustic  Modes 

Structural  characteristics  of  acoustic  modes  that  can  be 
generated  in  a  ramjet  are  of  considerable  importance  to  insta¬ 
bility  research,  particularly  with  regard  to  interactions  between 
the  acoustic  wave  and  combustion  and  fluid  dynamic  proces¬ 
ses  in  the  engine.  Mode  identification  is  possible  when  numer¬ 
ous  high-frequency-response  pressure  transducers  are  placed 
along  the  burner  (Tl  through  TIO,  Figure  I).  The  oscillatory 
pressure  data  from  the  10  transducers  are  recorded  on  multi¬ 
track  magnetic  tape.  Data  analysis  consists  of  a  series  of  tape 
playbacks  in  which  the  data  are  analyzed  for  frequency, 
amplitude,  and  phase.  The  results  of  a  large  number  of  tests  in 
four  different  combinations  of  inlet  length  and  combustor 
length  (Table  1)  have  shown  that: 

1 .  Four  different  acoustic  mode  structures  were  gener¬ 
ated. 

2.  More  than  one  acoustic  mode  structure  can  be  gener¬ 
ated  in  a  given  engine  configuration. 

3.  No  unique  relationship  exists  between  a  given  mode 
structure  and  fuel  type. 

The  results  of  the  acoustic  mode  structure  tests  are  sum¬ 
marized  in  Table  2.  The  third  mode  shape  in  that  table  (540  Hz) 
was  observed  in  the  majority  of  testing  to  date.  Its  characteris¬ 
tics  were  such  that  extensive  test  measurements  were  required 


Tabic  1 

Ramjet  Test  Configurations. 

Configuration 

no.a  Length,  in. 

Inlet 

Combustor 

1 

38.5 

30.1 

2 

38.5 

22.6 

3 

38  5 

15.1 

5 

28  5 

30.1 

•’Configuration  No  4  was  not  tested 


Tabla  2 

Observed  Acoustic  Modes 


MODE  SHAPE 


CONFIG  FREQ  FUEL 


NO  I  175  -  X 

NO  1  175  -  X  - 

NO  2  190  XXX 


(((((( 


(((((( 


s’x' 


NO  3  300 


NO  1  540 

NO  1  540 

NO.  5  540 


NO.  2  650  X  X 


combination  of  ethylene  and  JP-4  J 
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to  verify  its  structure.  The  test  measurements  were  also  com¬ 
pared  to  analytical  predictions  of  acoustic  mode  structure  that 
are  being  conducted  at  the  California  Institute  of  Technol¬ 
ogy.  11111 

The  analytically  predicted  pressure  amplitude  of  the  third 
mode  structure  (T able  2)  in  the  combustion  chamber  is  approx¬ 
imately  correct  for  a  half-wave  in  a  closed-closed  organ  pipe. 
In  contrast,  the  waveform  in  the  inlet  shows  the  strong  effects 
of  both  the  shock  and  the  mean  flow.  The  latter  causes  a 
reduction  of  the  wavelength,  a  result  of  the  Doppler  effect.  The 
shock  wave  acts  to  absorb  leftward-moving  waves,  a  process 
which  causes  the  waveform  to  move  upstream  during  a  cycle. 
This  appears  to  be  a  traveling  wave,  as  indicated  in  Figure  4, 
by  plotting  the  curve  at  t=0  and  again,  one-half  cycle  later. 
The  pattern  is  repeated  each  cycle  and  is  a  reflection  of  the  fact 
that  the  loss  of  acoustic  energy  at  the  shock  wave  introduces  a 
constant  phase  shift  between  leftward  and  rightward  moving 
waves  and  causes  the  reflected  wave  to  have  a  smaller 
amplitude  than  the  incident  wave.  The  wave  form  can  be 
characterized  as  a  nonstationary,  standing  wave. 

The  pressure  amplitude  distributions  in  Figure  5  show 
good  agreement  between  the  experimental  data  (data  points) 
and  the  analytical  treatment  (curve).  In  Figure  6,  the  relative 
phase  of  the  oscillations  (shown  at  a  later  time  during  one  cycle 
than  in  Figure  5)  are  shown  for  both  the  analytical  model 
(curve)  and  the  experiment  (data  points).  The  correlation  is 
gixxl;  the  exact  locations  of  the  phase  changes  are  different, 
but  all  the  major  features  coincide. 

In  general,  the  acoustic  mode  data  indicate  that  the 
mechanisms  that  cause  longitudinal  instabilities  in  the  research 
dump  combustor  have  a  relatively  broad  band  of  responsive¬ 
ness.  since  oscillations  have  been  observed  from  as  low  as  1 70 
Hz  to  as  high  as  650  Hz.  Further  details  concerning  results  of 
NWC  acoustic  mode  investigations  will  be  found  in  Refer¬ 
ences  19  and  20. 

Response  of  Inlet  Shock 

In  a  series  of  tests  to  determine  the  response  of  the  inlet 
shock  to  the  oscillating  diffuser  exit  pressure  (Pm,)  (540  Hz 
maximum),  the  oscillating  inlet  shock  position,  characterized 
by  the  diffuser  area  ratio  (diffuser  area  at  the  shock  position 
divided  by  area  of  the  throat.  ASH/A,),  was  measured  with 
Schlieren  photography  at  5.000  frames  per  second,  and  a 
time-average  shock  position  was  computed.  With  these  ex¬ 
perimental  data  the  time-average  shock  Mach  number  and  the 
time-average  Mach  number  in  the  diffuser  exit  (M,*)  were 
calculated.  By  comparing  these  time-average  experimental 
data  with  theoretical  data  based  on  one-dimensional  isen tropic 
theory  and  the  dimensions  of  the  ramjet  inlet,  the  subsonic 
diffuser  loss .  <* .  was  determined  as  a  function  of  Mw .  The  a  is 
defined  as  the  diffuser  exit  pressure  at  stagnation  divided  by 
the  diffuser  pressure  downstream  of  the  shock  at  stagnation. 


Figure  7 


Shock  displacement  amplitude  (peak  to-pcak)  versus 
diffuser  exit  pressure  amplitude. 
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By  applying  a  determined  for  the  steady-state  conditions 
to  the  oscillatory  conditions,  good  agreement  was  achieved 
between  the  theoretical  shock  displacement  amplitudes  (calcu¬ 
lated  from  2(PI)fc/PIx.)/l  —  (Pot/Pf^P  with  diffuser  exit  pres¬ 
sure  amplitude,  P„t,  and  mean  value  Pnt)  and  the  experimental 
shock  displacement  amplitudes  (ASH/A,)  (determined  from 
Schlieren  photography)  (Figure  7).-'1  This  agreement  indicates 
that  for  the  540-Hz  pressure  oscillation  the  shock  responded  in 
a  quasi-steady  fashion  to  the  pressure  fluctuation  amplitude. 

Summary 

The  results  of  the  long-term  program  to  determine  the 
nature  of  acoustic  combustion  instability  in  ramjets  have  been 
encouraging.  Predominant  frequencies  associated  w  ith  acous¬ 
tic  modes  of  the  ramjet  cavity  have  been  identified.  Good 
agreement  was  achieved  between  experiments  and  theory  for 
relative  pressure  amplitude  and  phase  distribution  in  the  ramjet 
cavity.  It  was  found  that  droplet  breakup/evaporative  phase  of 
liquid  fuels  is  not  a  necessary  condition  for  ramjet  combustion 
instability.  The  relationship  between  inlet  shock  and  diffuser 
exit  pressure  oscillations  and  the  chamber  acoustic  modes  and 
amplitudes  has  been  described.  Some  effective  applications  of 
one-dimensional  isentropic  theory  to  calculate  shock  motion 
amplitude  have  been  demonstrated.  Pursuit  of  the  implications 
of  these  and  future  findings  should  provide  solutions  to  ramjet 
combustion  instability  and  facilitate  deployment  of  tactical 
ramjets. 
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New  Concepts 
for  Thermal  Batteries 


by 

m.  H.  miles  and  A.  N.  Fletcher 
Naval  Weapons  Center 


Introduction 

The  first  thermal  batteries  were  developed  by  German 
scientists  during  World  Way  II  for  several  ordnance  applica¬ 
tions,  including  the  V- 1  and  V-2  missiles.1  Since  then,  thermal 
batteries  have  been  widely  used  in  military  applications,  in¬ 
cluding  various  missiles,  proximity  fuzes,  shells,  and  atomic 
weaponry.2  Maximum  power  from  a  thermal  battery  is  needed 
for  periods  ranging  from  a  few  seconds  to  60  minutes. 

As  electric  power  supplies  for  missile  weapon  systems, 
thermal  batteries  offer  many  advantages  over  alternate  energy 
sources.  These  advantages  include  long  shelf  life  (15-20 
years),  rapid  and  reliable  activation,  high  discharge  rates  ( 100 
mA/cm2),  high  power  densities,  and  the  ability  to  withstand 
high  shock  and  spin  forces.2  Many  of  these  advantages  result 
from  using  inert  salts  that  become  a  highly  conductive  electro¬ 
lyte  when  the  battery  is  activated  by  a  pyrotechnic  heat  source 
that  melts  the  salt. 12 

The  major  components  of  a  thermal  battery  cell  are  illus¬ 
trated  in  Figure  I .  The  anode  (negative  electrode)  and  cathode 
(positive  electrode)  in  the  galvanic  cell  are  separated  by  a  pad 
or  pellet  that  contains  the  depolarizer,  electrolyte  salt,  and 
binder.  Thermal  batteries  remain  inactive  and  inert  as  long  as 
the  salt  mixture  contained  between  the  anode  and  cathode 
remains  in  its  normal,  nonconducting,  inert,  solid  state.  The 
sudden  application  of  heat  activates  the  battery  by  melting  the 
salt,  converting  it  into  a  liquid  capable  of  conducting  electrical 
charges  (ions)  between  the  anode  and  the  cathode. 

Most  thermal  batteries  employed  in  weapon  systems  use 
a  lithium  chloride-potassium  chloride  (LiCI-KCI)  salt  mixture 
as  the  electrolyte  layer.  Calcium  (Ca)  metal  is  often  selected  as 


the  anode  material,  while  calcium  chromate  (CaCrO-i)  gener¬ 
ally  serves  as  the  active  cathode  material.  " 

The  use  of  a  LiCI-KCI  salt  mixture  has  several  disadvan¬ 
tages.  For  instance,  the  melting  point  of  the  mixture  limits  the 
battery's  activation  to  temperatures  above  350°  C.  Further¬ 
more.  the  Ca/LiCT-KCl/CaCrOi  system  exhibits  complex  side 
reactions  that  limit  the  operating  temperature  to  between  500 
and  600°  C  for  satisfactory  battery  performance. ' 4  These  high 
operating  temperatures  increase  the  quantity  of  pyrotechnic 
heating  material  and  insulation  required  for  cell  activation  and 
discharge.  Another  problem  with  LiCI-KCI  batteries  is  that 
CaCr04,  the  oxidizer  (depolarizer)  generally  used,  is  a  car¬ 
cinogen.'''  The  development  of  safety  regulations  regarding 
CaCr04  is  likely  to  curtail  its  availability  and  endanger  the 
production  of  thermal  batteries  that  use  it  as  an  oxidizer.  A 
promising  new  system  uses  molten  lithium  as  the  anode  and 
iron  disulfide  (FeS:)  as  the  oxidizer,  i.e..  Li/LiCI-KCl/FeS;. ' " 
This  system  is  free  of  carcinogenic  materials  and  is  superior  to 
the  Ca/CaCr04  system  in  power  and  energy  densities.  Fur¬ 
thermore,  this  new  thermal  battery  system  does  not  have  the 
performance-limiting  side  reactions  characteristic  of  the 
Ca/CaCrOj  battery  chemistries’ 

Molten  Nitrate  Systems 

The  novel  approach  of  using  certain  oxidizing  molten 
salts,  such  as  lithium  nitrate  (LiNO,),  in  thermal  batteries  has 
been  developed  at  the  Naval  Weapons  Center. 10  This  approach 
is  especially  attractive  because  of  the  low  melting  points  and 
nontoxic  nature  of  these  salts.  Mixtures  of  LiNO, -potassium 
nitrate  (KNO,)  salts  melt  at  temperatures  as  low  as  124°  C. 
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Figure  1 

Ma/or  components  of  a  simplified  thermal  battery  cell. 
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Figure  2 

Schematic  model  for  passivating  films  on  metals  The  term 
MXAV  represents  the  nonstoichiometric  film.  M ?•  represents 
the  excess  metal  ions  within  the  film,  and  X  represents  the 
anions  electrostatically  attracted  to  the  film-solution  interface 
Actual  films  will  not  necessarily  be  uniform  and  may  vary 
considerably  in  thickness. 
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Such  a  low  melting  temperature  can  shorten  a  battery’s  activa¬ 
tion  time  and  reduce  the  weight  of  heat  sources  and  insulation. 
The  LiNO,  salt  may  also  function  as  the  oxidizer,  thus 
eliminating  the  need  for  CaCrOj.  The  low  operating  tempera¬ 
tures  also  open  the  door  to  the  use  of  various  other  oxidizers, 
binders,  and  anode  materials  that  are  not  stable  at  higher  tem¬ 
peratures. 

Anode  Reactions 

Active  metals,  such  as  calcium,  that  are  used  as  the  ther¬ 
mal  battery  anode  are  not  normally  compatible  with  oxidizing 
molten  salts  such  as  LiNO,.  An  inert  protective  film,  however, 
limits  the  corrosive  chemical  reaction  between  the  anode  and 
the  molten  salt.  Regulation  of  this  passivating  film  is  essential 
to  good  anode  performance.  If  the  passivating  layer  becomes 
too  thick  and  hinders  the  passage  of  metal  ions  across  the  film 
and  into  the  solution ,  the  metal  ions  accumulate  within  the  film 
and  the  battery  performance  deteriorates.  For  the  calcium 
anode  in  molten  LiNO,,  aggressive  negative  ions,  such  as 
chloride  ions  (Cl  ),  must  be  added  to  the  molten  nitrate  to 
break  up  the  passivating  film  sufficiently  to  allow  the  thermal 
battery  cells  to  operate  efficiently. 11 

A  generalized  model  of  the  passivating  film  at  the  anode 
is  shown  in  Figure  2.  The  excess  metal  ions  (M7*)  trapped 
within  the  film  electrostatically  attract  negative  ions  (X  )  to 
the  film/solution  interface.  Passivating  negative  ions,  such  as 
nitrate  ions  (NO,  ),  peroxide  ions  (O,  ).  or  superoxide  ions 
(O,  ).  that  contain  oxygen  atoms  will  repair  breaks  in  a  pas¬ 
sivating  oxide  film.  Aggressive  negative  ions,  such  as  Cl  . 
hinder  the  film  repair  action  of  the  passivating  ions  through 
competitive  adsorption  at  the  film/solution  interface.12 

Theories  on  the  corrosion  of  metals  can  readily  be  tested 
in  molten  salt  electrolytes.  Corrosion  processes  are  much  more 
rapid  at  the  higher  temperatures  of  these  molten  salt  mixtures 
than  they  are  at  ordinary  ambient  conditions;  thus,  the  impor¬ 
tant  variables  can  be  readily  identified  and  studied.  Despite  the 
numerous  studies  of  corrosion,  there  is  no  general  recognition 
as  yet  of  a  specific  theory  or  model  that  explains  the  passivity 
and  corrosion  of  metals.  Generally  recognized,  however,  is 
the  concept  that  the  presence  of  chloride  ions  causes  pitting 
corrosion  of  metals  that  form  stable  oxide  films.  The  model 
presented  in  Figure  2  is  applicable  to  the  passivity  and  corro¬ 
sion  of  metals  under  a  variety  of  conditions,  including  the 
passivity  and  film  breakdown  of  active  metal  anodes  used  in 
thermal  battery  cells.12 

The  passivating  film  at  the  anode  hinders  the  passage  of 
metal  ions  into  the  solution;  thus  the  excess  of  positive  metal 
ions  trapped  within  the  film  causes  the  open  circuit  potential  of 
the  anode  to  become  more  positive.  With  the  breakdown  of  the 
passivating  film,  the  excess  metal  ions  escape  and  the  potential 
of  the  anode  shifts  to  more  negative  values.  Figure  3  shows  the 
breakdown  by  chloride  ions  of  the  passivating  film  on  calcium 
anodes  in  several  nitrate  melts.  Film  breakdown  was  deter¬ 
mined  by  the  sudden  change  in  the  open  circuit  potential  as 
measured  against  a  silver  ion-silver  (Ag*/Ag)  reference  elec¬ 
trode.  The  effect  of  chloride  ions  on  the  calcium  anode  poten¬ 
tial  is  much  greater  in  molten  LiNO,  than  in  either  molten 
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soldium  nitrate  (NaNO,)  or  molten  KNO,.  This  difference  is 
explained  by  considering  the  direct  reaction  of  molten  nitrates 
with  calcium  metal 

14  Ca  +  6  NO,  —  14  CaO|  +  2  N,|  +  N,Of  +  30 

to  produce  excess  oxide  ions  (O  )  in  the  melt .  These  oxide  ions 
are  innocuous  in  molten  UNO,  and  simply  precipitate  as  in¬ 
soluble  lithium  oxide  (Lip).  In  molten  NaNO,  and  KNO„ 
however,  oxides  are  more  soluble  and  form  the  more  reactive 
ionsO,  andO,  .Therefore,  according  to  the  model  presented 
in  Figure  2,  the  passivating  O,  and  0:  neagative  ions  formed 
at  the  calcium  surface  in  molten  NaNO,  and  KNO,  block  the 
desired  aggressive  action  of  the  added  Cl  ions.  As  a  result, 
smaller  shifts  are  observed  in  the  open  circuit  potential  of  the 
calcium  anode. 12  The  results  favor  the  use  of  UNO ,  rather  than 
NaNO, or  KNO ,  in  thermal  battery  cells  using  calcium  anodes. 

New  directions  for  research  on  the  anode  component  of 
thermal  battery  cells  are  suggested  by  the  passivating  film 
model  in  Figure  2.  The  nature  of  the  passivating  film  as  well  as 
the  rate  of  transport  of  metal  ions  across  the  film  would  be 
likely  to  vary  greatly  from  one  active  metal  to  another.  These 
considerations  suggest  the  study  of  other  active  metal  anodes, 
such  as  lithium  magnesium,  strontium,  scandium,  and  yt¬ 
trium,  in  molten  nitrates.  Because  of  its  much  smaller  size,  the 
lithium  ion  (Li  *)  has  a  faster  transport  rate  across  the  passivat¬ 
ing  film  than  do  larger  metal  ions  such  as  Ca**.  The  low 
melting  point  of  lithium  metal  ( 1 88°  C),  however,  requires  that 
this  metal  be  used  in  the  form  of  higher  melting  compounds  or 
alloys.  Preliminary  studies  have  shown  that  the  use  of  lithium 
composites  as  anodes  in  molten  nitrates  yields  high  current 
densities  with  very  little  polarization. 

Cathode  Reactions 

The  second  major  component  of  a  thermal  battery  cell  is 
the  cathode.  Oxidizing  molten  salts  such  as  nitrates  can  func¬ 
tion  both  as  the  electrolyte  and  the  active  cathode  material .  The 
performance  of  the  molten  nitrate  cathode  is  determined  prin¬ 
cipally  by  the  nature  of  the  positive  metal  ions  present.  Ex¬ 
perimentally,  the  potential  of  the  LiNO ,  cathode  is  about  a  full 
volt  better  than  that  of  the  KNO,  cathode.14  The  performance 
of  the  molten  nitrate  cathode  varies  inversely  with  the  metal 
ion  radius  and  decreases  in  the  order  LiNO,  >  NaNO,  > 
KNO,.  A  model  for  the  electrochemical  reduction  of  molten 
nitrates  is  presented  in  Figure  4.  The  stronger  positive  electric 
field  surrounding  the  smaller  metal  ion  aids  both  the  slow 
transfer  of  the  electron  from  the  electrode  to  the  nitrate  ion ,  and 
the  rupture  of  the  N-O  bond. I4  |:' 

The  cathodic  reaction  in  molten  LiNO,  cells 

LiNO,  +  2  Li*  +  2e  — »  LiNO:  +  Lip  i 

produces  insoluble  Lip  as  a  product.  Severe  cell  polarization 
occurs  when  Lip  is  formed  at  the  electrode  surface  faster  than 
it  can  be  removed,  thus  blocking  the  reacting  LiNO,  from 
reaching  the  electrode  surface  and  slowing  the  rate  of  the 
cathodic  reaction.  The  polarization  of  the  cathode  becomes  the 
limiting  factor  in  obtaining  useful  cell  voltages  at  current  den¬ 
sities  above  50  mA/cm2. 


Figure  3 

Effects  of  addition  of  10-mole-percent  chloride  salts  on  the 
open-circuit  potentials  of  the  calcium  anode  in  molten  KNO 
NaN03.  and  UN03  at  350 °  C.  Chloride  additions  were  made 
at  time  zero 
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Cathode  Performance  Improvement 

Adding  small  amounts  of  silver  nitrate  (AgNO,)  to  mol¬ 
ten  nitrates  greatly  improves  the  performance  of  the  cathode."’ 
Figure  5  shows  that  adding  0. 1 -mole-percent  AgNO,  to  mol¬ 
ten  LiNO,  at  350°  C  permits  the  cathode  to  sustain  a  current 
density  of  100  m  A/cm2  without  the  onset  of  the  severe  passiva¬ 
tion  observed  in  pure  LiNO,.  The  deposition  of  silver  particles 
on  the  electrode  surface  increases  its  effective  area,  thereby 
minimizing  the  blocking  action  caused  by  the  insoluble  Li  p 
formed.  In  the  presence  of  I  0-mole-peicent  AgNO,.  the  di¬ 
rect  reduction  of  silver  ions  ( Ag ' )  yields  stable  potentials  close 
to  that  of  the  Ag‘/Ag  reference  electrode  (Figure  5). 2 

Long-term  discharge  experiments  indicate  that  a  recy¬ 
cling  of  silver  ions  takes  place  by  reactions  such  as 

Ag*  +  e  -»  Ag 

2  Ag  +  LiNO,  +  2Li*  -*  2Ag*  +  LiN02  +  Lip|. 

This  corrosive  oxidation  of  the  deposited  silver  metal  by  the 
molten  LiNO,  recycles  the  Ag  *  ions.  This  reaction,  like  many 
others  in  molten  nitrates,  depends  on  the  nature  of  the  positive 
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Figure  4 

Postulated  reaction  steps  for  the  electrochemical  reduction  of 
molten  nitrates  showing  the  involvements  of  the  positive 
metal  i on 
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metal  ion  present  and  is  much  slower  in  molten  NaNO,  or 
KNO,  than  in  molten  LiNO,. 

Various  other  metal  ion  additives  have  also  been  tested  in 
molten  nitrates.  Of  those  tested,  only  lead  ( Pb  * ' )  and  cad¬ 
mium  (Cd*  *)  ions  showed  any  evidence  for  a  cathodic  reac¬ 
tion;  however,  in  contrast  to  the  reaction  of  silver  ions,  there  is 
no  evidence  for  a  new  metallic  phase  ever  forming  on  the 
electrode  surface.  Naval  Weapons  Center  studies  suggest  that 
reduction  to  the  free  metal  never  occurs  for  Pb  *  •  or  Cd’  *  ions 
in  molten  nitrates;  instead  the  NO,  ion  in  Pb(NO,)  *  or 
Cd(NO,)*  complexes  is  the  actual  substance  reduced;  i.e., 

M(NO,)*  +  2e  — >  MO  +  NO,  . 

Unlike  the  cathodic  reaction  of  silver  ions,  electrons  rather 
than  metal  ions  transfer  the  charge  across  the  electrode/molten 
nitrate  interface.  The  strong  electric  field  of  the  divalent  Pb  *  * 
or  Cd**  ion  in  these  complexes  lowers  the  free  energy  of 
activation  for  the  reduction  of  the  nitrate  ion.  Thus  a  new 
reduction  wave  at  a  more  positive  potential  is  observed.  These 
studies  substantiate  the  major  role  of  the  metal  ion  on  the 
electrode  reduction  of  molten  nitrates  (Figure  4). 
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Figure  5 

Effect  of  silver  nitrate  addition  on  the  performance  of  the 
cathode  during  constant  current  studies  at  100  mA:cm?  in 
molten  UN01  at  350°  C 


Cell  Discharge 

Results  of  the  separate  fundamental  studies  of  the  anode 
and  cathode  in  molten  nitrates  can  be  used  to  design  thermal 
battery  cells  that  yield  the  optimum  performance.17  A  thermal 
battery  system  using  these  results  is  represented  by  the  cell 

Ca/LiCI  .LiNO , /LiNO, .AgNO,/Ni 

with  the  LiNO,  incorporated  into  a  fiberglass  binder.  LiCI 
layered  near  the  anode,  and  AgNO,  layered  near  the  cathode. 
The  fiberglass  binder  also  acts  as  a  separator  and  confines  the 
LiCI  and  AgNO,  additives  to  the  anode  and  cathode  regions, 
respectively,  where  they  are  most  effective.  A  typical  dis¬ 
charge  of  this  ceil  is  shown  in  Figure  6.  At  a  discharge  rate  of 
75  m  A/cm2 ,  the  cell  potential  remains  above  2  volts  for  over  60 
seconds.  This  cell  was  designed  to  supply  a  high  current  den¬ 
sity  for  a  short  time.  Other  thermal  battery  cells  that  use  molten 
nitrates  can  be  designed  with  different  properties,  such  as  dis¬ 
charge  times  exceeding  60  minutes  at  10  mA/cm!. i;-16 
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Figure  1 


Introduction 


Modern  military  explosives 


RDX  MMX 


The  next  generation  of  naval  ordnance  now  in  the  plan¬ 
ning,  research,  and  development  cycle  requires  new  energetic 
materials  to  fulfill  new  ordnance  performance  and  survivabil¬ 
ity  requirements.  To  meet  these  requirements,  the  Naval 
Weapons  Center  supports  an  active,  multifaceted,  in-house 
research  effort  in  energetic  materials  synthesis  and  evaluation. 
Programs  range  from  research  projects  of  theoretical  interest  in 
the  design  and  synthesis  of  new  explosive  molecules,  to 
applied  projects  on  formulation  and  evaluation  of  new  explo¬ 
sive  and  propellant  systems.  The  program  is  guided  by  the 
requirements  defined  by  the  systems  design  engineers  and 
includes  new  high-density  high-energy  explosives  for 
volume-limited  systems  and  insensitive  high  explosives  for 
applications  that  call  for  new  explosives  that  are  insensitive  to 
shock  and  extreme  temperatures.  All  new  systems  need  to 
minimize  hazards  associated  with  manufacturing,  processing, 
loading,  deployment,  and  demilitarization.  The  research  goal 
at  N  WC  is  to  combine  these  diverse  needs  and  develop  the  new 
generation  of  energetic  materials  required  by  future  weapon 
systems. 


Modern  Military  Explosives 

Most  modem  military  explosives  are  composed  of  a  bin¬ 
der  material,  a  solid  crystalline  explosive  compound,  and  var¬ 
ious  additives.  The  binder,  generally  a  polymer,  has  two  func¬ 
tions:  it  holds  the  composite  mixture  together,  and  it  imparts 
mechanical  properties  to  the  mixture.  Additives  are  used  to 
adjust  the  chemical,  physical,  and  explosive  properties  of  the 
resultant  explosive  formulation.  The  solid  crystalline  explo¬ 
sive  compounds  most  often  used  in  high-performance  explo¬ 
sive  formulations  are  2.4,6-trinitrotoluene  (TNT),  1.3,5- 
trinitro- 1 ,3.5-triazacydohexane  (RDX),  1 ,3,5,7-tetranitro- 
1 ,3.5,7-tetraazacyclooctane  (HMX),  and  1 ,3,5-triamino- 
2,4, 6-trinitrobenzene  (TATB).  (See  Figure  I  for  the  chemical 
compositions  of  these  compounds  and  Table  1  for  the  physical . 
chemical,  and  explosive  properties  of  these  compounds.) 


Tabla  1 

Physical.  Chemical,  and  Explosive  Properties 
ot  i  NT.  TATB,  RDX.  and  t<MX 

Melting  Detonation  Detonation  Impact 


Compound 

Density  point 

(gccl  (°C) 

velocity 

(ms) 

pressure 

(kbar) 

sensitivity 
(Hw  cm) 

TNT 

1  65 

81 

6.930 

210 

75 

TATB 

1  94 

325 

7.760 

291 

350s 

RDX 

1  81 

205 

8.700 

338 

28 

HMX 

1  90 

285 

9.110 

390 

20 

*  Estimated 
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If  the  next  generation  of  energetic  materials  is  to  be  use¬ 
ful,  the  materials  should  have  the  density  and  energy  of  RDX 
and  HMX  and  the  reduced  sensitivity  to  an  inadvertent  initia¬ 
tion  of  TNT  and  TATB .  These  extreme  requirements  may  be 
mutually  exclusive,  but  reasonable  tradeoffs  of  performance 
and  hazard  properties  can  and  must  be  made  to  develop  the 
required  materials.  Research  is  providing  new  methodology 
and  new  materials  for  future  use  in  energetic  materials  applica¬ 
tions,  which  include  polynitroaromadc,  insensitive,  and  new 
polynitramine  explosives. 

Polynitroaromatic  Explosives 

Hexanitrobenzene  (HNB)  belongs  to  a  class  of  zero  hyd¬ 
rogen  explosives  called  nitrocarbons  and  is  composed  only  of 
nitro  groups  attached  to  carbon  (C,(NO.O,).  After  a  more  than 
40-year  effort  by  free  world  scientists,  NWC  research 
chemists  succeeded  in  1978  in  synthesizing  HNB  from  pen- 
tanitroaniline  (PNA)  as  shown  in  Figure  2.'  HNB  is  the  most 
energetic  organic  explosive  known.  The  key  to  the  HNB 
synthesis  was  the  development  of  a  new,  generally  applicable 
procedure  for  the  oxidation  of  weakly  basic  amines  to  nitro 
functionality.  This  new  procedure  uses  in  situ  generated 

(peroxydisulfuric  acid  as  the  oxidant.  HNB  was  only  the  third 

nitrocarbon  explosive  reported  in  the  literature.  The  synthesis 
procedure  was  used  to  prepare  many  new  polynitro-substitutcd 
j  explosives  and  was  recently  adapted  for  the  preparation  of  a 

J  new  zero  hydrogen  explosive,  decanitrobiphenyl  (DNBP)  as 

I  shown  in  Figure  3.- '  DNBP  joins  the  select  list  of  zero  hydro¬ 

gen  nitrocarbon  explosives,  along  with  HNB,  tetranit- 
romethane  (TNM),  and  hexanitrocthane  (HNE).4 '  DNBP  was 
prepared  from  4,4'-diamino-2,  2', 3,  3',5,5',6,  6'-octanitro- 
!  biphenyl  (ONBP)  by  oxidation  with  peroxydisulfuric  acid, 

j  ,  DNBP  is  chemically  similar  to  but  more  reactive  than  HNB. 

[  Present  observations  suggest  that  the  thermal  stabilities  of  the 

!  known  nitrocarbons  decrease  in  the  order  of  tetranitromethane 

>  hexanitroethane  >  hexanitrobenzene  >  decanitrobiphenyl. 
Research  efforts  are  continuing  in  this  area  of  explosives 
synthesis. 

Insensitive  High  Explosives 

Energetic  materials  research  at  NWC  also  involves  the 
synthesis  of  new  insensitive  high  explosives  having  both  exp¬ 
losive  properties  equal  to  or  greater  than  RDX  or  HMX  and 
sensitivity  properties  similar  to  those  of  TATB  in  resistance  to 
thermal,  shock,  and  impact  initiation.  The  amino  group  has 
been  identified  as  a  mediator  of  explosive  sensitivity.  It  also 
serves  as  an  activating  group  for  electrophilic  aromatic  sub¬ 
stitution  reactions.  The  presence  of  an  amino  group  on  the 
aromatic  nucleus  thus  serves  a  dual  purpose  of  allowing  the 
simultaneous  introduction  of  nitro  groups  in  the  aromatic  nuc¬ 
leus  and  imparting  the  desired  sensitivity  properties.  This  dual 
nature  is  illustrated  by  the  synthesis  of  DNBP  as  shown  in 
Figure  3.  ONBP  is  prepared  by  nitration  of  the  tetranitro  pre¬ 
cursor,  4,4'-diamino-2,2',6,6'-tetranitrobiphenyl  (TNBP). 


Figure  2 

Synthesis  of  HNB 


Figure  3 

Synthesis  ot  decanitrobiphenyl. 
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This  nitration  would  not  proceed  without  the  activation  and 
directive  influence  of  the  4.4'  amino  groups.  Oxidation  of 
ONBP  amino  groups  is  accomplished  using  peroxydisulfuric 
acid.  Comparison  of  the  ONBP  versus  the  DNBP  impact  sen¬ 
sitivities  shows  the  sensitivity  mediation  effects  of  the  4.4' 
amino  substituents. 

The  amino  group  with  its  dual  nature  was  used  advan¬ 
tageously  to  prepare  a  large  number  of  previously  inaccessible 
aromatic  explosives.  TNT  has  served  as  a  very  convenient 
precursor  to  these  new  explosives.  TNT.  which  is  available  in 
vast  quantities  at  low  cost,  can  be  selectively  reduced  to  yield 
two  isomeric  dinitrotoluidines.  Ammonium  hydroxide 
catalyzed  reduction  of  TNT  by  HS  in  dioxane  solvent  results 
in  para-nitro  group  reduction  giving  4-amino- 
2 ,6-dinitrotoluene  (p-DNT).  The  ortho  isomer.  2- 
amino-4,6-dinitrotoluene  (o-DNT).  is  obtained  by  reduction 
of  TNT  with  iron  metal  in  acetic  acid.  The  remaining  meta 
isomer  (m-DNT)  could  not  be  obtained  directly  from  TNT.  but 
was  obtained  from  meta  toluic  acid  (MTA)  in  several  steps. 
The  three  dinitrotoluidines  can  be  readily  nitrated  in  mixed 
acid  to  give  the  desired  tetranitrotoluidines  as  summarized  in 
Figure  4.  The  three  isomeric  tetranitrotoluidines  have  better 
calculated  explosive  properties  than  TNT.  The  densities  are 


slightly  higher  than  TNT  but  less  than  expected  when  calcu¬ 
lated  by  the  Holden  method  (1.76  g/cc).a  The  tetrani¬ 
trotoluidines  are  slightly  more  sensitive  to  impact  than  TNT. 
but  considerably  less  sensitive  than  RDX  or  HMX  as  shown  in 
Table  2. 


Table  2. 

Comparative  Properties  of  Tetranitrotoluidines. 

Detonation  Detonation  Impact 
Density  velocity  pressure  sensitivity 


Compound 

(g  cc) 

(ms) 

(kbar) 

(H50  cm) 

TNT 

1  65 

6,900 

210 

70 

Ortho 

1  72 

7,700a 

269a 

36 

Meta 

T  73 

7.700a 

269a 

37 

Para 

1  72 

7.70 0a 

269a 

47 

HMX 

1.90 

9.100 

390 

25 

d  Calculated  ' 

Figure  4 

Synthesis  of  the  isomeric  tetranitrotoluidines 
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The  three  isomeric  tetranitrotoiuidines  can  be  used  as 
effective  precursors  of  other  explosives.  Oxidation  using 
peroxydisulfuric  acid  gives  pentanitrotoluene  (PNT)  in  high 
yields.  The  overall  conversion  of  TNT  to  PNT  is  80%  when 
prepared  from  para-tetranitrotoluidine,  making  this  excellent 
explosive  readily  available.  PNT  has  a  measured  density  of 
1.76  g/cc  (calculated  =  1.82  g/cc),  has  very  good  calculated 
explosive  properties  (detonation  velocity  =  8,200  m/s,  deto¬ 
nation  pressure  =  312  kbar) ,  and  exhibits  sensitivity  similar  to 
RDX  (22  cm  versus  25  cm). 

When  the  mixed  acid  nitrations  of  the  dinitrotoluidines 
are  carried  out  at  reduced  temperatures  and/or  acid  strength, 
pemitrated  intermediates  can  be  isolated,  as  shown  in  Figure 
5,  for  the  para  isomer.  These  nitramine  intermediates  are 
thermally  unstable  and  slowly  evolve  nitrogen  oxides  if  al¬ 
lowed  to  stand  at  room  temperature.  However,  if  the  nitramine 
is  suspended  in  concentrated  sulfuric  acid  and  then  treated  with 
anisole,  the  nitramine  is  cleaved,  giving  high  yields  of  the 
corresponding  tetranitrotoluidine.  Attempts  to  purify  these  in¬ 
termediate  nitramines  by  recrystallization  led  to  the  discovery 
that  they  undergo  facile  thermal  rearrangement  to  yield  ben- 
zenediazooxides  as  shown  in  Figure  5  for  PNToL.  This  facile 
nitramine  rearrangement  is  characteristic  for  polynitro  aroma¬ 
tic  nitramines  and  has  made  possible  the  preparation  of  many 
new  extremely  sensitive  and  powerful  polynitroben- 
zenediazooxides.  (Diazooxides  are  known  as  primary  explo¬ 
sives:  DINOL,  4,6-dinitrobenzene-2-diazo-l-oxide  was  used 
as  a  primary.)  These  new  explosives  may  be  useful  as  lead 
azide  replacements  in  new  ordnance. 

Para-tetranitrotoluidine  has  been  found  to  undergo  un¬ 
ique  chemistry  in  its  own  right.  When  the  mixed  acid  nitration 
of  the  precursor  dinitrotoluidine  (p-DNT)  is  carried  out  at 
elevated  temperatures,  a  remarkable  transformation  occurs — 
not  only  is  the  ring  pemitrated,  but  the  methyl  substituent  is 
also  removed  and  replaced  by  a  fifth  nitro  group  to  give  50  to 
55%  yield  of  pentanitroaniline  as  shown  in  Figure  6.  This 
unexpected  transformation  involves  oxidative  demethylation 
of  the  aromatic  ring  followed  by  normal  mixed  acid  nitration  of 
the  vacated  para  position.  The  importance  of  this  discovery  is 
that  PNA  serves  as  a  convenient  precursor  to  two  very  impor¬ 
tant  explosives:  HNB  and  TATTJ.  Oxidation  of  PNA  with 
peroxydisulfuric  acid  gives  high  yields  of  HNB.  If  PNA  is 
dissolved  in  methylene  chloride  and  then  treated  with  am¬ 
monia,  high  yields  of  TATB  are  obtained.  Readily  available 
TNT  becomes  a  precursor  not  only  for  the  most  powerful 
organic  explosive  yet  synthesized  (HNB)  but  also  to  the  ar¬ 
chetypal  insensitive  high  explosive  TATB .  The  rich  chemistry 
of  TNT  is  shown  in  Figure  7. 


Figure  S 

Nitramine  rearrangement  to  diazooxide  and  cleavage  by 
anisole  to  the  free  amine. 
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Conversion  of  p-DNT  to  PNA 
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Synthesis  of  Nitramine  Explosives 

As  noted,  earlier,  the  cyclic  aliphatic  nitramine  explo¬ 
sives,  RDX  and  HMX,  are  the  most  energetic  materials  now 
used  in  high-performance  high-energy  applications  of  naval 
ordnance  and  propellant  systems  The  performance  of  new 
explosives  is  measured  against  these  benchmark  materials. 
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Figure  7 

Explosives  derived  from  TNT.  HE  is  high  explosive  and  IHE  is 
insensitive  high  explosive 


Pf  N  f  ANI  T  HO  AN'L  INI 
(Hi) 


Tt  TRANITROTOIUIDINES 
<HE) 


HI  HIGH  EXPIOSIVI 

INSENSITIVITY  HIGH  EXPLOSIVE 


Empirical  calculations,  practice,  and  experience  indicate  that 
RDX  and  HMX  are  not  at  the  limit  of  energy  content  achieva¬ 
ble  from  this  class  of  high-explosive  compounds.  Calculations 
show  that  increasing  ring  size  in  these  compounds  may  lead  to 
increased  density.8  RDX  and  HMX  can  be  considered  the 
second  and  third  members  of  a  homologous  series  of  nitramine 
explosives.  Other  members  of  this  series  are  unknown  but  in 
principle  can  be  prepared  simply  by  adding  an  additional 
methylene-nitramine  unit  to  the  ring  skeleton  as  shown  in 
Figure  8. 

Constraining  the  molecular  structure  further  by  additional 
bridging  to  yield  polycyclic  molecules  also  may  induce  in¬ 
creased  energy ,  density,  and  stability.  (Density  will  increase  in 
the  structural  alteration:  acyclic  <  monocyclic  <  bicyclic  < 
polycyclic  <  condensed  polycyclic  (caged)  molecule.)  Typi¬ 
cal  molecules  of  this  type  are  shown  in  Figure  9  and  are  rep¬ 
resentative  of  the  target  compounds  of  the  nitramine  explosive 
synthesis  program  now  being  carried  out  at  NWC. 

The  synthesis  of  macrocyclic  and  polycyclic  nitramine 
explosives  depends  on  the  development  of  new  synthetic 
techniques  and  methodologies.  Recent  work  has  concentrated 
on  developing  this  new  technology  and  applying  it  to  the 
synthesis  of  new  nitramine  explosives.  A  new  method  for  the 


Figure  8  Figure  9 

Macrocyclic  nitramine  explosives.  Cyclic  nitramine  explosive  target  molecules. 
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synthesis  of  1 ,3-dinitro- 1 ,3-diazacycloalkanes  was  de¬ 
veloped.9  The  discovery  that  alpha  omega  diamines  can  be 
trapped  in  situ  with  formaldehyde  to  form  the  corresponding 
dinitroso  compound,  when  treated  with  nitrous  acid,  was 
pivotal.  Nitrolysis  of  the  dinitroso  intermediates  proceed  in 
high  yield  to  give  the  corresponding  nitramine,  when  reacted 
with  100%  nitric  acid  as  shown  in  Figure  10.  Reaction  of  an 
alpha  omega  diamine  with  formaldehyde  generates  a  1,3- 
diazacycloalkane  (aminal).  The  aminal  is  trapped  in  situ  by 
adding  nitrous  acid  yielding  a  1 ,3-dinitroso- 1 ,3- 
diazacycloalkane.  Nitrolysis  proceeds  with  excellent  yields 
using  100%  nitric  acid  to  give  the  desired  nitramine.  This 
technique  has  been  applied  to  the  synthesis  of  several  interest¬ 
ing  polycyclic  polynitramine  compounds  such  as  2,4,6, 8- 
tetranitro-2 ,4 ,6 ,8-tetranitrospiro [5 . 5Jundecane  (TNSU) . 
TNSU  was  prepared  from  tetrakis(aminomethyl)methane 
(TAM)  as  shown  in  Figure  1 1 .  TNSU  is  more  energetic  than 
TNT  yet  has  approximately  the  same  impact  sensitivity  (65 
cm). 

Members  of  the  decalin  family  of  tetraazanitramines 
were  prepared  using  similar  methodology:  trans- 1 ,3,5,7- 
tetranitro-l,3,5,7-tetraazadecalin  was  prepared  from  1 ,2,3,4- 
tetraamino  butane,  which  was  derived  in  several  steps  from 
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Figure  12 

Synthesis  of  as-  and 
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diethyltartrate.  The  isomeric  c/'j-l,3,5,7-tetranitro-l  ,3,5,7- 
tetraazadecalin  was  prepared  in  a  similar  manner  (see  Figure 
12).'° 

Substitution  of  glyoxal  for  formaldehyde,  which  is  an 
adaptation  of  this  new  procedure,  was  used  to  prepare  another 
member  of  this  family  of  energetic  materials,  irons- 1 ,4,5,8- 
tetranitro-l,4,5,8-tetraazadecalin."  This  material  was  pre¬ 
pared  by  condensation  of  glyoxal  with  ethylene  diamine  fol¬ 
lowed  by  nitrosation  to  give  the  tetranitroso  derivative.  Nit- 
rolysis  of  the  nitroso  substituents  was  accomplished  using 
100 %  nitric  acid  as  shown  in  Figure  13.  The  three  decalins  are 
all  well-behaved  solids  with  good  thermal  stability.  Their 
physical  and  explosive  properties  are  listed  in  Table  3. 


Table  3. 

Properties  of  the  Tetranitrotetraaza  dec  aim  Series 


Compound 

Melting 

point 

CC) 

Density 

(g/cc) 

Detonation  Impact 
velocity  *  sensitivity 
(m/s)  (H^,  cm) 

Trans-2,4,5,8 

234 

1.80 

821 

40 

CIS- 1, 3.5,7 

238 

1  78 

8.21 

Trans-1.3,5,7 

252 

1.75 

821 

a  Calculated 


The  methodology  described  above  has  been  used  to  pre¬ 
pare  a  new  explosive  compound  with  properties  that  exceed 
those  of  RDX  and  HMX.  The  synthetic  strategy  employed 
required  the  replacement  of  ethylene  diamine  with  a  suitable 
energetic  1,2-diamine,  4,5-diaminofurazan  (DAF)  for  the 
condensation  reaction  with  glyoxal.  Modified  literature 
methods  were  used  for  the  diamine  preparation. i:  "  14  As 
shown  in  Figure  14,  when  the  diamine  is  condensed  with 
glyoxal.  a  good  yield  of  the  desired  bis-furazanodecalin  (BFD) 
parent  compound  is  obtained .  The  nitration  of  the  tetramine  by 
nitric  acid  in  trifluroacetic  acid  (TFA)  proceeds  as  expected, 
and  CL- 15  ( 1 ,4,5,8-tetranitro-l  ,4,5,8-tetraazadifurazano 
(3.4-c)(3,4-h)decalin)  is  prepared  in  good  yield.15  A  com¬ 
parison  of  the  properties  for  CL- 15  and  HMX  is  shown  in 
Table  4. 
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Properties  of  CL-1 5  and  HMX 
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Density 

velocity 
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Compound 
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CL- 15 
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9.57a 
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1.90 
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^Calculated  .,s 


Figure  14 

Synthesis  of  CL- 15. 
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In  the  early  studies  of  CL-  IS  the  compound  was  noted  to 
be  thermally  unstable.  Samples  of  the  compound  left  at  room 
temperature  soon  evolved  red  gas;  samples  left  for  prolonged 
time  periods  decomposed  to  intractable  material.  The  heat  of 
formation  of  CL- 15  is  +185  kcal/mole  as  determined  by  the 
standard  back  calculation  technique  from  the  heat  of  combus¬ 
tion.  Using  the  heat  of  formation  and  the  measured  density,  the 
performance  parameters  were  calculated  by  the  Kamlet 
method.'6  CL-15  exceeds  HMX  in  explosive  performance. 
While  CL- 1 5  does  not  have  the  required  thermal  properties ,  its 
successful  synthesis  demonstrates  the  possibility  of  exceeding 
HMX  performance.  Preparation  is  expected  of  other  explo¬ 
sives  having  this  energy  content  and  the  requisite  stability 
properties. 
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Introduction 

Many  physical  processes  within  the  outermost  layers  of 
atoms  in  solids  differ  from  those  deeper  within  the  atoms'  bulk 
material.  Chemical,  electronic,  optical,  mechanical,  and  other 
processes  are  changed  in  detail  or  effect  by  the  asymmetrical 
surface  environment  of  the  atom  or  by  interaction  with  the 
external  ambient.  For  example,  the  ( 1 10)  plane  of  a  GaA.x 
crystal  contains  Ga-As  bonds  lying  within  the  plane  when  the 
plane  is  within  the  bulk.  If  the  crystal  is  cleaved  in  vacuum 
along  the  (110)  plane,  the  atoms  rearrange  themselves  so  that 
the  Ga-As  bond  tilts  out  of  the  plane  by  about  28  degrees. 
Electronic  changes  occur  along  with  this  "reconstruction"  of 
the  (110)  surface  plane.  The  absence  of  neighboring  atoms  on 
the  vacuum  side  frees  bonding  or  valence  orbitals  on  each  atom 
so  that  the  orbitals  "dangle"  in  the  vacuum.  These  surface 
dangling-bond  states  are  eigenstates  of  the  reconstructed  two- 
dimensional  surface  lattice,  but  they  are  not  eigenstates  of  the 
bulk  crystal.  Because  they  are  not  part  of  the  bulk  band  struc¬ 
ture.  they  may  act  to  scatter  or  trap  electrons  from  Bloch 
eigenstates  of  the  three-dimensional  cry  stal.  Charges  trapped 
in  such  surface  or  interface  states  are  sources  of  electrostatic 
potential  that  affects  the  motion  of  charge  carriers  w  ithin  hun¬ 
dreds  of  lattice  spacings  from  the  surface.  In  this  way.  the 
action  of  the  outermost  surface  layer  of  an  atom  may  influence 
charge  transport  or  optical  properties  in  a  much  larger  region  of 
the  solid. 

The  study  of  the  outermost  lay  ers  of  atoms  is  the  domain 
of  surface  science.  Many  techniques  are  exploited  in  the  study 
of  surface  or  interface  layers.  Such  layers  may  be  recon¬ 
structed  as  a  consequence  not  only  of  vacuum  cleavage,  but 
also  of  adsorption,  oxidation,  metal  film  deposition,  heteroepi- 
taxial  growth,  or  other  surface  modifications.  Most  surface- 
science  techniques  involve  stimulus-response  measurements, 
using  a  probe  beam  to  stimulate  the  surface  layer.  Usually  ,  the 
surface  response  is  observed  by  the  emitted  or  scattered  parti¬ 
cles.  Such  experimental  techniques  may  be  categorized  by  the 
probe-beam  particles  and  the  detected  particles,  either  of 
which  may  be  ions  (or  atoms),  electrons,  or  photons  Ideally  . 
the  probe  beam  should  excite  only  one  atomic  layer  with  a 
single-quantum  excitation,  without  causing  significant  dam 
age  or  change  in  the  surface.  Ion  probe  beams  penetrate  least 
deeply  into  the  solid,  exciting  the  fewest  layers  but  causing  the 
greatest  damage  to  the  surface.  Photon  probe  beams  penetrate 
much  more  deeply .  exciting  many  layers  of  atoms  but  stimulat¬ 
ing  mainly  electronic,  single-quantum  excitations  and  causing 
little  or  no  atomic  disorder  in  the  surface.  Electron  probe 
beams  penetrate  quite  deeply,  causing  some  atomic  disorder 
while  inducing  a  variety  of  electronic  excitations  by  inelastic 
scattering.  No  probe  beam  approaches  the  ideal  of  single¬ 
layer,  single-quantum  excitation  without  irreversible  distur¬ 
bance  to  the  surface. 

Ideally  ,  detected  particles  should  carry  information  about 
an  undisturbed  single  layer  of  surface  atoms  or,  even  more 
specifically  ,  about  a  particular  atomic  site  in  the  surface  layer. 
Generally,  the  emitted  low-energy  ions  represent  the  outer¬ 
most  surface  layer,  while  the  electrons  or  photons  emitted 


represent  electronic  excitations.  Again,  none  of  the  detected 
particles  ideally  represents  specific  excitations  of  the  surface 
atomic  layer  alone,  but  combinations  of  the  probe-beam  parti 
cles  and  the  detected  particles  provide  valuable  information 
about  the  surface  layer  of  atoms. 

Just  as  point  defects  in  the  bulk  crystal  may  control  or 
dominate  the  electrical  or  optical  properties  of  cry  stals,  point 
defects  on  the  crystal's  surface  may  strongly  influence  the 
electrical  and  optical  properties  of  the  surface.  For  example, 
plots  of  optical  absorption  versus  thickness  usually  show  a 
nonzero  intercept,  representing  an  absorption  associated  with 
the  surfaces.  Usually  this  absorption  is  ascribed  to  imperfec¬ 
tions.  including  chemical  impurities,  in  the  surf  ace  layers  that 
arc  not  distributed  throughout  the  bulk.  Electronic  or  vibra¬ 
tional  excitations  of  surface  imperfections  are  observed  in  sur¬ 
face  optical  absorption  measurements.  However,  detailed 
theory  of  surface  optical  absorption  is  in  a  primitive  state  be¬ 
cause  of  the  difficulty  off  measuring  and  calculating  the  elec¬ 
tronic  states  associated  with  surface  point  defects.  Only  re¬ 
cently  .  surface-science  measurement  techniques  that  are  sensi¬ 
tive  to  such  localized  states  in  surface  layers  have  become 
available.  Two  of  these  techniques  are  surface-extended, 
x-ray -absorption  fine  structure  (SEXAFS)  and  photon- 
stimulated  ion  desorption  ( PSID* ).  Both  of  these  techniques 
provide  sensitivity  to  the  outermost  surface  layer,  and  PSID 
also  provides  sensitivity  to  isolated  point  defects. 

SEXAFS  Techniques 

The  use  of  ordinary  extended  x-ray  absorption  fine  struc¬ 
ture  (EXAFS)  techniques  does  not  provide  surface  sensitivity 
because  of  the  deep  penetration  of  x-rays  into  the  crystal. 
However,  in  SEXAFS  techniques  surface  sensitivity  is  pro¬ 
vided  by  detecting  the  photoelectron  yield  (PEY).  instead  of 
the  transmitted  x-rays  resulting  from  soft  x-ray  excitation.  By 
detection  of  PEY.  the  x-ray  absorption  within  the  photoelec¬ 
tron  escape  depth  (4- 10  A)  is  monitored  as  a  function  of  photon 
energy.  Careful  analysis  of  SEXAFS  data  provides  near¬ 
neighbor  interatomic  distances  around  the  surface-excitation 
site  For  example,  if  the  spectral  region  above  the  As  ( M ,).  or 
2s.  absorption  edge  is  studied  on  a  GaAs  (1 10)  surface,  the 
Ga-As  bond  length  and  possibly  the  As- As  separation  may  be 
obtained  on  the  reconstructed  surface.  Such  a  study  may  be 
made  as  the  surface  is  exposed  to  oxygen,  sub- monolayer 
metal  coverage,  elevated  temperature,  ionizing  radiation,  or 
other  stimuli  that  may  influence  the  interatomic  distances  in 
the  surface  region. 

In  principle.  SEXAFS  measurements  could  be  used  to 
probe  the  env  ironments  of  isolated  atoms  on  surfaces  ( such  as 
Au  atoms  on  GaAs).  or  specific  groups  of  photoelectrons  (such 
as  angle-  or  energy-resolved  photoelectrons)  could  be  detected 
in  place  of  the  total  PEY  In  practice,  the  SEXAFS  signal  is  t<x> 
weak  to  permit  observation  with  the  reduced  signal  strength 

*  PSID"  in  sometimes  referred  to  as  '  PSD"  in  analog) 
with  electron  stimulated  desorption  fl-.SD) 
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Figure  1 

PSID  results, m  a  co- evaporated  NbsSn  thin  film  alter  a  5 00°C 
bake  m  vacuum  Part  A  shows  the  time-ol-thght  mass 
spectrum  ot  photoemitted  positive  ions.  Both  molecular  and 
atomic  surface  mipurties  may  be  observed  Part  B  shows  a 
representative  excitation  spectrum  By  observing  that  P 
emission  shows  a  strong  threshold  at  the  Si(L  j  excitation 
energy  three  levels  ol  knowledge  result.  ( 1 1  analysis:  Si  and  F 
are  present  on  the  Nb  ,Sn  surface.  (2)  chemistry  F  bonds  to  Si 
mpurity  atoms  on  the  Nb  ,Sn  surface.  (31  physics  the  spectral 
shape  of  the  near-edge  structure  may  be  analyzed  to  study 
me  electronic  structure  at  the  Si-F  surface-impurity  site  on 
Nb  iSn  i  see  reference  41 


these  measurements  would  require.  SEXAFS  is  used  in  the 
study  of  sites  of  majority  speeies  on  the  surface,  but  not  to 
study  the  sites  of  minority  species.  To  study  the  structure  and 
electronic  states  at  sites  of  minority  species  on  the  surfaces, 
PSID  may  he  used.  The  remainder  of  this  paper  will  deal  with 
PSID  and  its  applications. 

PSID  Techniques 

Exploiting  the  single-quantum  excitation  characteristic  of 
the  photon-probe  Seam  and  the  single-layer  sensitivity  of  emit¬ 
ted  low  -energy  ions.  PSID  techniques  facilitate  the  study  of  a 
variety  of  structural  and  electronic  features  of  isolated  surface 
sites.  In  PSID  measurements,  low  -energy  ions  that  are  emitted 
or  "desorbed"  from  the  surface  are  identified  by  charge-to- 
mass  ratio  using  a  time-of-flight  or  other  mass  analyzer. 1  The 
excitation-energy  spectrum  for  each  ionic  speeies  emitted  is 
obtained  by  using  a  scanning  monochromator  in  the  photon- 
probe  beam.  Empirically .  ii  has  been  found  that  the  probability 
of  PSID  is  very  low.  below  10-15  eV.  and  that  photon- 
excitation  thresholds  for  ion  emission  coincide  with  the  core- 
i  xcihiiion  energies  of  surface-site  atoms.-  Thus.  PSID  has 
been  shown  to  result  from  a  single-quantum  excitation  of  an 
atomic-core  electron. 

Because  the  kinetic  energy  of  the  emitted  ion  is  so  low 
'/.<  10  eVi.  it  is  concluded  that  the  wins  issue  from  the  outer¬ 
most  surface  layer.  The  emission  of  higher-energy  ions  after 
energy  loss  hy  collisions  is  ruled  out  by  time-resolved  results 
and  by  the  emitted  ion  kinetic-energy  spectrum.  Time- 
resolved  results  show  that  any  delay  time  for  an  energy -loss 
process  must  be  less  than  10  ''  sec  and  must  be  independent  of 
excitation  energy  above  the  threshold  energy  .  Emitted-ion 
kinetic-energy  spectra  are  generally  peaked  and  do  not  show 
the  marks  of  a  collisional  energy -loss  process.  Thus  PSID 
measurements  detect  ions  from  the  outermost  layer  of  surface 
atoms  as  stimulated  by  a  single-quantum  atomic-core  excita¬ 
tion. 

Isolated  sites  in  the  surface  layer  may  he  studied  by  PSID 
in  two  ways.  Isolated-surface  impurity  atoms  may  he  iden¬ 
tified  directly  among  desorbing  speeies,  or  indirectly  as  sur¬ 
face  sites  for  other  desorbing  species  Both  methods  of  iden¬ 
tification  are  needed  for  a  complete  study  of  all  the  chemical 
impurities  on  surfaces  because  of  the  wide  variation  in  the 
desorption  probability  for  various  atoms.  Singly  bonded 
atoms,  such  as  Group  I  or  Group  VII  atoms,  are  readily  de¬ 
sorbed  while  multiply-bonded  atoms  are  seldom  desorbed. 
Multiply-bonded  atoms  arc  studied  as  PSID  sites.  With  two 
complementary  methods,  many  surface  impurities  can  he 
studied  using  PSID  techniques.'1 

Surface  Chemical  Analysis  of 
Superconducting  Thin  Films 

Electronic  devices  using  superconductors  have  de¬ 
monstrated  extremely  fast  switching  times  (<  10  11  s).  Super¬ 
conducting  microwave  cavities  have  extremely  low  hiss  ( 1/Q 
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<  10  g).  Superconducting  magnetometers  are  capable  of  de¬ 
tecting  magnetic-field  changes  <10  12  T.  With  such 
capabilities,  it  seems  remarkable  that  superconducting  devices 
are  not  more  widely  used  in  both  commercial  and  military 
applications.  Problems  impeding  the  broader  application  of 
superconducting  devices  involve  both  the  difficult  materials 
technology  and  the  costly  cryogenic  technology.  If  reliable 
superconducting  devices  can  be  fabricated  from  high-T 
superconductors  such  as  NbiGe  or  Nb,Sn,  both  cost  and  per¬ 
formance  factors  will  improve;  materials-technology  de¬ 
velopment  is  the  keystone. 

Traditionally,  superconducting  materials  have  been  con¬ 
sidered  sensitive  to  metallurgical  treatment,  which  affects 
crystal  structure,  but  insensitive  to  the  trace  chemical  im¬ 
purities  that  dominate  semiconductor  properties.  On  the  other 
hand,  the  Josophson-junction  device  derives  its  dramatic  per 
formance  from  the  interaction  of  superconducting  charge  ear¬ 
ners  (Cooper  pairs)  with  a  thin  insulating  barrier.  Tunnel  cur¬ 
rents  through  the  barrier  depend  on  barrier  parameters,  mainly 
thickness  and  height.  The  barrier  height,  in  turn,  depends  on 
the  chemical  composition  and  purity  of  the  extremely  thin 
tunnel  barrier,  which  is  usually  <50  A  thick.  Thus,  in 
Joxephson-junction  or  similar  devices,  minor  impurities  on 
superconductor  surfaces  may  play  an  important  role  in  the 
performance  of  the  device. 

The  first  PSID  study  of  chemical  impurities  on  Nb  and 
Nb  ,Sn  surfaces  clearly  illustrated  the  complex  surface  chemis¬ 
try  involved.4  In  all.  25  impurities  were  identified  by  PSID. 
while  only  carbon  and  oxygen  were  identifiable  by  Auger 
electron  spectroscopy.  Hydrogen  and  oxygen  are  expected  to 
be  of  greater  significance  in  tunnel-junction  performance, 
especially  if  the  oxidized  amorphous-silicon  (a-Si)  barrier  is 
used.  Further  studies  are  being  carried  out  on  the  roles  of  these 
impurities,  both  in  the  growth  of  high-quality  Nb,Ge  films  and 
in  the  formation  of  tunneling  barriers  of  oxidized  a-Si  on 
high-T  superconducting  films  (Figure  I). 


Figure  2 

Atomic  Sites  on  the  (1 10)  surface  of  Group  lll-V  compound 
semiconductors.  The  surface  lattice  shown  (bold lines)  is  the 
(1)0)  plane  of  the  truncated  crystal  tot  zinc-blende  structure) 
Group  III  atom  sites  are  represented  by  solid  circles  ana 
Group  V  atom  sites  on  the  surface  are  represented  by  open 
circles  In  the  "reconstruction  "  of  the  clean  surface,  all  Group 
III  surface  atoms  sink  slightly  into  the  crystal,  and  all  Group  V 
surface  atoms  rise  slightly  out  of  the  crystal  This 
reconstruction  results  in  the  tilling  of  the  Group  lll-V  surface 
bond  out  ot  the  surface  plane  by  about  28  degrees  Each 
surface  atom  has  one  free  bond  dangling  out  of  the  crystal 
because  of  the  absence  ot  the  next  atomic  layer  Electron:  : 
shifts  cause  the  dangling  bond  o>  the  Group  III  atom  to  be 
unoccupied  (above  the  Fermi  energy),  its  electron  shifting  tc 
till  tht:  dangling  bond  ot  the  Group.  V  atom,  which  lies 
energetically  deep  m  the  valence  band 


Hydrogen  Interaction  With  (1 10)  Surfaces 
of  Group  III-V  Compound 
Semiconductors 

As  mentioned  in  the  i .production,  electronic  surface 
states  on  semiconductors  may  trap  or  scatter  charge  carriers 
even  though  the  charge  earners  are  in  Bloch  eigenstates  of  the 
three-dimensional  crystal.  For  this  reason,  the  control  of  sur¬ 
face  electronic  states  is  important  in  electronic  devices,  and 
increasingly  important  as  the  device  speed  is  increased  or  the 
device  size  is  reduced.  One  of  the  simplest  chemical  agents  for 
surface-state  control  is  hydrogen,  which  forms  a  single,  simple 
covalent  bond  with  many  atomic  species.  Although  the  room- 
temperature  reactivity  of  hydrogen  with  Group  lll-V  semicon¬ 
ductor  surfaces  is  low,  the  high  sensitivity  of  PSID  to  hydro¬ 
gen  facilitates  the  study  of  very  small  amounts  of  surface 
hydrogen  (Figure  2). 
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Figure  3 


H'  PSIDexcitatr  ,n  spectra  of  cleaved ( 1 10>  GaAs  Hydrogen 
a,  ss:  '/n'i/  'ii  G<iAs  during  crystal  growth  is  bonded  to  Ga 
sort, ret-  atoms  on  freshly  cleaved  surfaces  The  excitation 
rrsOolifs  shown  near  20  eV  correspond  to  Ga  tM.:  ..i  or  3d 
>  -XC  tat"  ’ti  at  the  sudds  e  Ga-H  Site  Fxpos'ng  the  surface  to 
H  O  .at  sir  t'r,  u/i«  cs  ar  add'ts trial  H '  desorption  thrt -sholdat 
33  eV  irresrondiiig  to  OlL  •  '  or  2s  ex- labor  ot  the  surface 
OH  d\ is  irpate  (see  rv’ereriste  b> 


The  study  of  the  hydrogen  interaction  with  surface  dan¬ 
gling  bonds  begins  with  the  study  of  freshly  cleaved,  clean 
semiconductor  surfaces.  Group  lll-V  semiconductors  cleave 
along  (110)  cry  stal  planes.  Such  clean  surfaces  are  prepared  in 
situ  in  ultrahigh  vacuum,  where  surfaces  contaminate  slowly 
over  a  period  of  hours  by  interaction  with  ambient  gas 
molecules  PSII)  on  such  freshly-cleaved,  clean  surfaces 
show  s  little  or  no  ion  y  ield  because  the  multiply  -bonded  Group 
lll-V  elements  are  seldom  desorbed  An  exception  was  found 
for  GaAs  crystals  grown  in  a  hydrogen  environment.  The 
freshly  cleaved  surfaces  of  these  crystals  showed  considerable 
H  desorption  when  the  Ga(M,,f.  or  fd.  core  states  were 
excited'  (f  igure  4). 

The  interaction  of  hydrogen  with  InP  ( 1 1(1)  surfaces  was 
studied  using  PSII)  and  total  PHY  A  low  pressure  stream  of 
excited,  molecular  hydrogen  (H  ‘I  was  flushed  across  the 
freshly  cleaved  surface  to  an  exposure  level  of  KIM..  PSII) 
and  total  PHY  excitation  spectra  were  then  studied  above  the 
P<  K  i  and  Im  1.)  excitation  thresholds."  The  total  PHY  spectrum 
was  nearly  unchanged  by  the  H ■*  exposure,  while  the  PSID 
results  showed  H  desorption  spectra  with  thresholds  at  P(kl 
and  at  the  three  (ml.)  edges,  f  sing  the  total  Pf-.Y  spectra  near 
these  edges  as  a  measure  of  the  photoabsorption  in  the  surface 
region,  it  was  deduced  that  the  population  of  In-H  surface 
bonds  and  P-H  surface  bonds  is  in  the  same  order  of  mag¬ 
nitude  .  Prev ions  studies  of  the  GaAs  (110)  surface  indicated  a 
similar  result  Thus,  the  reactivity  of  H  *  with  Group  III  and 
Group  V  elements  on  ( 1 10)  surfaces  appears  to  be  approxi¬ 
mately  equal  (Figure  4). 

Ihis  result  suggests  a  model  for  hydrogen-surface  in¬ 
teraction  involving  the  dangling-bond  states.  The  H.  molecule 
has  two  electrons,  and  two  more  electrons  are  found  in  the  P 
dangling  bond,  facilitating  the  formation  o (adjacent  In-H  and 
P-H  bonds.  At  the  same  time,  the  establishment  of  similar 
electron  densities  around  the  In  and  P  surface  atoms  by  the 
adsorption  of  hy  drogen  removes  the  cause  of  the  bond  tilt,  and 
the  In-P  bond  is  expected  to  return  to  the  (110)  plane  as  in  the 
bulk  crystal  (Figure  5). 

The  hy  drogen  dosage  used  in  these  experiments  (  UP  L) 
w  as  far  less  than  needed  to  completely  cover  the  surface  w  ith 
hydrogen.  Although  the  hydrogen  coverage  is  unknown,  it  is 
expected  to  be  less  than  0.01  monolayer,  possibly  much  less. 
Published  research  on  the  interaction  of  oxygen  with  GaAs 
under  similar  circumstances  led  to  conclusion  that  a  dose  of 
more  than  10"'  L  of  O.*  would  be  necessary  to  completely 
cover  the  surface.  Our  results  on  hydrogen  chemisorbed  on 
InP,  GaP.  and  GaAs  surfaces  probably  represent  isolated  hyd¬ 
rogen  surface  sites  or  small  covered  areas.  To  test  the  hydrided 
surface  for  transport  properties,  we  must  learn  to  produce 
stable,  complete  coverage  of  the  surface.  If  such  coverage  of 
the  surface  can  be  done,  surface  trapping  and  scattering  of 
charge  carriers  at  hydrided  Group  Ill- V  semiconductor  sur¬ 
faces  may  be  expected  to  be  reduced  because  the  surface  dangl¬ 
ing  bonds  are  used  in  covalent  bonds  similar  to  the  bulk  bonds . 

+  l.  =  tjngmuir.  a  unil  ol  exposure  equivalent  to  10  " 
lore  seconds. 
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Figure  4 

Total  PEY  and  PSID  excitation  spectra  compared  at  the 
P(K)  and  ln(Lt  edges  on  hydrogen-dosed  InP  (110)  (see 
reference  6) 
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Figure  5 

Model  tor  hydrogen  surface  sites  on  InP  The  left-hand  In-P 
bond  is  shown  tilted  as  appropriate  to  a  clean  InP  (1 10) 
surface  with  surface  dangling  bonds  On  the  right,  hydrogen 
is  shown  bonded  to  ad/acent  In  and  P  surface  atoms  Both 
dangling  bonds  are  now  occupied  and  the  electronic 
symmetry  restores  the  bulk-crystal  positions  of  the  atoms  (see 
reference  6) 
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The  ffifchelson 
Laboratory  Synchroton 
Radiation  Project 

by 

Victor  Rehn 
Naval  Weapons  Center 


The  Michelson  Laboratory  Synchrotron  Radiation  Pro¬ 
ject  ( MLSRP)  originated  at  the  Nava]  Weapons  Center(NWC) 
in  1972  to  exploit  synchrotron  radiation  (SR)  in  the  study  of 
electronic  materials  that  have  potential  for  Navy  application. 
In  the  decade  since  1972,  many  types  of  SR  experiments  have 
been  conducted  involving  nearly  two  dozen  NWC  scientists. 
Each  of  the  remarkable  properties  of  SR  has  been  exploited: 
the  high-intensity,  highly  polarized,  low-noise  continuum 
spectrum  (visible  to  x-ray),  the  ultraclean,  ultrahigh  vacuum 
environment,  the  subnanosecond  radiation  pulses  with  a 
1.28-MHz  repetition  rate,  and  the  highly  collimated  small 
beam.  MLSRP  continues  under  NWC  and  Office  of  Naval 
Research  (ONR)  funding,  with  current  emphasis  on  surface 
studies  using  photon-stimulated  ion  desorption  described  in 
another  article,  “Surface  Science  With  Photon-Stimulated  Ion 
Desorption."  in  this  issue. 

The  project  had  its  origin  in  August  1972  at  a  Gordon 
Conference  on  high-energy  spectroscopy,  where  interest  was 
in  band-structure  analysis  for  wide-gap  semiconductors  and 
insulators.  Professor  F.  C.  Brown  of  the  University  oflllinois 
mentioned  SPEAR,  the  newly  built  Stanford  positron -electron 
acceleration  ring,  as  a  potentially  wonderful  source  of  SR  for 
vacuum-ultraviolet  (VUV)  spectroscopy.  Later  Professor  W. 
E.  Spicer  of  Stanford  University  agreed  that  a  collaboration  in 
the  construction  of  experimental  facilities  and  subsequent  ex¬ 
periments  on  VUV  reflectance  and  modulation  spectroscopies 
would  be  welcomed.  In  subsequent  meetings  with  manage¬ 
ment  personnel  at  NWC  and  Stanford,  the  basis  for  collabora¬ 
tion  and  funding  proposals  was  formulated.  In  November 
1972,  a  proposal  was  submitted  by  NWC  requesting  funding 
for  a  joint  NWC-ONR  project  in  VUV  SR  research  in  elec¬ 
tronic  materials. 


In  early  February  1973.  ONR  funding  was  received.  The 
project  began  in  earnest.  Professor  S.  Doniach  of  Stanford  and 
Dr.  Spicer  formed  similar  collaborations  w  ith  Dr.  Brown,  then 
at  the  Xerox  Palo  Alto  Research  Center,  and  with  Professor  J. 
Baldeschwieler  of  the  California  Institute  of  Technology.  In 
June  1973,  the  Stanford  SR  project  (SSRP).  was  funded  by  the 
National  Science  Foundation  (NSF)  on  the  basis  that  NSF 
would  fund  the  building  construction,  the  control  and  safety 
systems,  and  the  modifications  necessary  in  the  SPEAR  ring. 
Each  of  the  collaborating  groups  was  responsible  for  one  beam 
line  and  was  to  be  granted  priority  for  initial  use  of  that  beam 
line  after  official  opening  of  the  experimental  program.  NWC 
was  responsible  for  the  VUV  beam  line  (now  called  the 
8-degree  branch  line  or  beam  line  1-2)  and  mirror  design  and 
procurement  for  the  soft  x-ray  beam  line  built  by  the  Xerox 
group  (Figure  I).  In  November  1973.  P.  Eisenberger  of  Bell 
Laboratories  joined  the  collaboration  to  build  an  extended 
x-ray  absorption  fine  structure  (EXAFS)  station.  By 
November  1974,  all  five  experiment  stations  were  function¬ 
ing,  and  an  exciting  variety  of  SR  experiments  at  SSRP*  was 
beginning. 

Among  the  MLSRP  accomplishments  of  the  past  eight 
years  was  the  measurement  of  the  VUV  transmission  and  re¬ 
flection  spectra  of  lithium  yttrium  fluoride.1  LiYF4.  ami  the 


•“SSRP”  in  now  called  “Stanford  Sy  nchrotron  - 
Radiation  Laboratory  (SSRLt.”  and  is  a  national  facility  for 
research  usinj*  synehrolrtvn  radiation,  bunded  by  IV  pari  mem 
ol  l.neruy.  SSRI.  is  available  to  all  qualified  users,  and  has 
been  used  by  more  than  7()0  scientists. 
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Figure  1. 


Beam  lines  of  the  original  SSRP  This  complex  of  five 
branching  beam  lines  is  now  referred  to  as  beam  line  1 
Although  MLSRP  constructed  only  the  VUV  beam  line, 
currently  MLSRP  is  a  "user"  of  all  SSRL  beam  lines 
appropriate  to  the  experimental  requirements 
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VUV  reflectance  and  scattering  spectra  of  various  SiC  mate¬ 
rials.2  The  former  experiment  showed  Li  YF4  to  have  a  funda¬ 
mental  absorption  edge  at  10.5  eV  ( I 185  At,  making  it  com¬ 
petitive  with  LiF  and  MgF.  among  UV  transmitting  materials. 
The  SiC  results  helped  to  establish  its  remarkable  optical  prop¬ 
erties  and  led  to  its  grow  ing  use  for  mirrors  or  mirror  sub¬ 
strates.  Other  achievements  were  the  stud>  of  fluorescence 
excitation  and  emission  spectra  of  matrix-isolated  small 
molecules  by  time-resolved  spectroscopies,’1  and  the  study  of 
VUV  and  soft  x-ray  scattering  of  real  mirror  surfaces.* 

The  highlight  was  the  discovery  of  photon-stimulated  ion 
desorption  (PSID)  in  December  1 D7H M.  L.  Knotek  of  San- 
dia  Laboratories  had  presented  the  Knotek-Feibclman  model 
for  electron-stimulated  ion  desorption  from  fully  valent  ionic 
solids.  In  this  model ,  the  first  step  is  creation  of  a  core  hole  in  a 
surface  atom  by  inelastic  electron  scattering.  Knotek  and  I 
were  convinced  that  the  subsequent  events  leading  to  ion  de¬ 
sorption  were  independent  of  the  mechanism  of  creation  of  the 
core  hole.  I  suggested  that  this  idea  could  be  tested  using  SR. 
and  that  the  NWC  apparatus  for  angle- resolved  scattering  to¬ 
gether  with  the  time-resolved  spectroscopy  electronic  system 
could  be  adapted  for  the  test.  I  pointed  out  the  advantages  of 
time-of-flight  mass  spectroscopy  for  noise  reduction  and  for 
identification  of  emerging  ionic  species.  Seventy-two  hours  of 
beam  time  in  December  1978  were  "bootlegged"  for  the  test. 
The  predictions  were  confirmed.  H  ' .  OH  * ,  and  O  ions  were 
observed  to  desorb  from  a  TiO,  crystal  surface,  after  reducing, 
ion  etching,  and  exposing  the  surface  of  H  :()  vapor.  "’  Since  that 
time  PSID  has  been  applied  to  the  study  of  semiconductors, 
insulators,  metals,  superconductors,  molecular  solids,  and  ad¬ 
sorbed  layers.  A  wide  variety  of  scientific  phenomena  has 
been  investigated,  including  surface  chemical  analysis,  sur¬ 
face  bonding,  photoioni/.ation  processes,  adsorbate-surface 
chemical  reactions,  surface  dangling-bond  charging,  and 
photolysis. 

MLSRP  has  also  produced  important  results  of  a  less 
tangible  nature.  For  example,  attention  was  drawn  to  the 
potentialities  and  technical  problems  of  free-electron  lasers  for 
military  applications  and  to  certain  special  materials,  such  as 
SiC  for  optical  components  for  high-energy  lasers.  Closer  in¬ 
teractions  were  fostered  between  NWC  and  university  or  in¬ 
dustrial  laboratories,  both  here  and  abroad.  These  interactions 
include  interdisciplinary  ones  that  grew  from  the  heterogene¬ 
ous  mix  of  specialists  who  use  SR  in  their  research.  This 
symbiotic  laboratory  environment  may  tum  out  to  have  as 
great  a  benefit  to  the  U.S.  technological  community  as  will  the 
scientific  discoveries  made  in  the  SR  laboratory. 
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The  mathematics  of 
Simple  Random  Walks 

by 

Daniel  T.  Gillespie 
Naval  Weapons  Center 


Introduction 

Many  real-world  processes  contain  some  element  of  ran¬ 
domness.  When  this  random  element  is  so  strong  that  the 
course  of  the  process  is  not  uniquely  determined  by  knowable 
parameters,  the  process  is  said  to  be  nondeterministic .  The 
best  known  examples  of  a  nondeterministic  process  are  games 
of  chance,  such  as  craps  or  roulette.  Some  examples  of  a  more 
pressing  concern  to  the  Navy  are  the  depletion  of  weapons 
stores  on  a  combat  ship  and  the  evolution  of  a  complex  war 
game.  Nondeterministic  processes  also  abound  in  the  physical 
sciences.  For  example,  since  we  cannot  track  all  molecular 
positions  and  velocities  in  a  gas  of  many  thermally  colliding 
molecules,  it  follows  that  the  evolution  of  a  gas-phase  chemi¬ 
cal  reaction  is  basically  a  nondeterministic  process. 

In  a  deterministic  process,  the  state  S(t)  of  the  system  at 
any  time  t  suffices  to  uniquely  determine  the  state  SU  +  tit)  at 


an  infinitesimally  later  time  t  +  at.  This  fact  allows  construc¬ 
tion  of  a  differential  equation  for  the  state  function  5 .  and  if  the 
state  of  the  system  at  some  initial  time  is  specified,  the  corres¬ 
ponding  solution  of  that  differential  equation  unambiguously 
gives  the  state  at  any  later  time.  But  if  the  process  is  nondeter¬ 
ministic  —  e.g.,  if  a  knowledge  of  SU)  permits  only  a  prob¬ 
abilistic  prediction  to  be  made  of  SU  +  dt)  —  this  traditional 
analytical  approach  clearly  docs  not  work. 

The  analysis  of  nondeterministic  processes  has  been  the 
focus  of  a  small  but  persistent  investigative  effort  in  the  inde¬ 
pendent  research  program  at  the  Naval  Weapons  Center.  One 
of  the  main  goals  has  been  simply  to  team  and  develop  basic- 
mathematical  approaches  and  techniques.  In  that  somewhat 
tutorial  spirit,  we  shall  review  here  some  of  the  mathematical 
ways  of  analyzing  one  simple  kind  of  nondeterministic  pro¬ 
cess.  Although  our  presentation  contains  much  that  is  well 
known,  it  also  incorporates  many  clarifications  and  new  re¬ 
sults  developed  in  our  research  program. 
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A  Simple  Random  Walk  Model 

We  consider  a  system  whose  possible  states  can  be  put 
into  a  one-to-one  correspondence  with  the  non-negative  inte¬ 
gers  x.  Furthermore,  we  suppose  that  the  time  evolution  of  the 
system  is  governed  by  two  given  functions,  a  ,(.v)  and  a  Ur), 
according  to  the  following  scheme: 

a  .(.r  )dt  =  the  probability  that  the  system,  if  in  state  x  at 
time  /,  will  jump  to  state  ,v±l  in  the  next 
infinitesimal  time  interval  (t.t+dt). 
U=0.l.2.  ...)  (I) 

We  shall  assume  that  a  (0) =0,  so  that  the  system  stays  off  the 
negative  v-axis. 

One  example  of  such  a  system  is  a  mixture  of  chemically 
reacting  molecules  in  which  the  populations  of  all  but  one 
reactant  are  held  constant;  in  that  case,  x  gives  the  current 
number  of  molecules  of  the  unconstrained  species,  and<i  ,(.v) 
and  a  (.r)  are  the  probabilities  per  unit  time  of  chemical  reac¬ 
tions  which,  respectively,  create  and  destroy  one  of  those 
molecules.  Another  example  is  a  bound  cluster  of  .v  molecules 
floating  in  a  vapor  of  single  molecules;  in  that  case,  a  ,(.v)  and 
a  (v)  are  the  probabilities  per  unit  time  for  the  cluster  to  con¬ 
dense  and  evaporate  a  molecule. 

Mathematically,  our  model  is  an  example  of  a  ‘  contin¬ 
uous-time.  discrete-state  Markovian  process."  The  term 
“Markovian"  conveys  the  fact  that  the  random  behavior  of  the 
system  is  independent  of  its  past;  this  is  true  because  the  prob¬ 
abilities  in  equation  ( I )  depend  at  most  on  the  current  state  x 
and  not  on  how  the  system  got  to  that  state.  F.quation  ( I ) 
implies  that  the  state  of  the  system  moves  about  on  the.c -axis  in 
sudden,  irregular  jumps  from  one  integer  to  another.  It  is 
customary  to  refer  to  such  behavior  as  a  "random  walk,"  and 
to  the  system  itself  as  the  "walker."  Our  particular  random 
walk  is  further  constrained  to  be  a  "single- variable,  single- 
step"  random  walk,  since  the  state  of  the  walker  is  specified  by 
a  single  variable  that  can  change  by  only  one  unit  at  a  time. 

The  Master  Equation 

One  way  to  mathematically  describe  our  random  walk 
model  is  to  focus  on  the  function 

Pix.t)  =  the  probability  that  the  system  will  be  in  the 

state  x  at  time  t .  (2) 

We  can  obtain  a  time-evolution  equation  for  P  by  first  express¬ 
ing  P(x,t  +dt)  as  the  sum  of  the  probabilities  of  the  three  possi¬ 
ble  routes  to  being  in  state  x  at  time  t+dt: 

Pix.t  +dt)  =  P(x  -  I  ,t)  •  a  ,(.r  -  I )  dt  +  P(.r  +  1 .1 ) 

•a  ( x  +  l)dr  +  P(x,r)  •  (I  -a  .(xidr-a  (,vW/|. 

Here,  the  first  term  on  the  right  is  the  probability  that  the 
system  is  in  state  x  - 1  at  time  t  and  then  jumps  up  one  state  in 
the  next  dr,  the  second  term  is  the  probability  that  the  system  is 
in  state  x  + 1  at  time  /  and  then  jumps  down  one  state  in  the  next 
dt;  the  third  term  is  the  probability  that  the  system  is  already  in 


state  r  at  time  t  and  then  does  not  jump  at  all  in  the  next  dt.  A 
simple  algebraic  rearrangement  leads  to  the  following  time- 
evolution  equation  for  the  function  P: 

i>  P(x,t)l b  t  =  P(.v -  I  ,t)a  Ax- I )  +  P(.r  +  I  j)a  (x  +  I )  - 

P(x,t)  («,(.v)+<i  tv)].  (3) 

This  differential-difference  equation  is  called  the  master  equa¬ 
tion.  Although  it  is  easy  to  derive,  it  is  usually  extremely 
difficult  to  solve. 

If  /  is  any  function  of  the  state  x.  then  the  average  or 
expected  value  of/  at  time  i  is 

(/'(.»)>,  =  i/(.v» Plx.t). 

\  0 

The  n'h  moment  of  P  is  defined  to  be  (.t’ ),.  Of  particular 
interest  are  (x  >,,  the  average  of  v  at  time  t,  and  |  <.r  >,  - 
tv  >,“ ] 1  the  standard  deviation  of  .v  at  time  t;  the  latter  is  a 
measure  ot  the  expected  difference  between  an  observed  value 
of.v  at  timer  and  the  average  (x  >,.  Multiplying  equation  (3)  hv 
.<»  and  then  summing  over  v.  we  obtain  the  following  time- 
evolution  equation  for  the  moments  of  P: 

d  {x"),ldt  =  (a  Ax)  It  v  +  I )"  -  .»■"]),  - 

{a  (vl|.v”-u  -  lrj),  (4) 

For  n  =0  this  equation  show  s  that  the  time  evolution  specified 
by  equation  (3)  is  consistent  with  the  normalization  require¬ 
ment.  2. /*(.«- ./)  =  I.  For /i  =  I  equation  (4)  gives 

d (x  ),ldl  =  (a  ,(x)— a  t.v )  >, .  ( 5 1 

For  n  a:  I  we  have  a  problem  w  ith  equation  (4)  if  either  a  << )  or 
a  Ax )  is  nonlinear  in  x.  In  that  case  the  right  side  will  contain 
tv'"),  for  some  m >n ,  thereby  rendering  the  set  of  equations 
infinitely  openended.  Unfortunately,  no  known  general  pro¬ 
cedure  exists  for  solving  such  a  set  of  equations. 

If  the  functions  a  .(.» )  and  a  (  t )  are  every  w  here  positiv  e 
(except  that  a  (<)>=(>)  and  satisfy  a  Ax)la  (< )— «•()  as  i— ►  *. 
then  it  can  be  shown  that  P(i.r)— *Pj.r)  as  t— **,  where  the 
function  P  „  is  independent  of  both  /  and  the  initial  state  of  the 
walker,  (This  of  course  does  not  imply  that  the  w  alker  eventu¬ 
ally  stops  walking.)  Equation  (3)  evidently  requires  that  P, 
satisfy 

0  =  Pax  -  I  )a  ,t»  -  I )  +  Pax  +  I  )a  U  +  I ) 

-  P*U)l«  ,(r)  +  a  tv)]. 

A  simple  transposition  of  terms  reveals  that  the  quantity 
P.tv)o  tv)  -  P.tv-l)tf.tv-l)  must  be  a  constant  for  allx. 
and  taking  v  =0  shows  that  this  constant  must  be  zero.  Thus. 

P  Ax)  =  (u  ,(r  -  I  da  (.v))Pj(r-l)  U  =  l,2,...)  (6) 

This  is  a  simple  recursion  relation  for  P„  iteration  of  which 

gives 

P Ax)  —  P AO)  11  |<i  An  -  I  )la  (n)|.  (v  =  1.2 _ )  (7) 

n  I 

The  quantity  P  jO)  may  be  obtained  from  equation  (7)  and  the 
normalization  condition.  i,P  J.v)=  I . 
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Figure  1 

The  function  Pjx)  for  a  bistable  system .-'s 
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The  extrema  of  the  function  Px  may  be  found  (to  within 
one  unit  of  t)  by  setting  its  •’derivative"  [PJ.* )-PJjc-  1 )]/ 1 
equal  to  zero  and  then  solving  for  *.  Using  equation(6),  we  can 
easily  show  that  the  extrema  of/*  .are  those  points.*  for  which 

a  Ax- I )-«  (.*)  =  0.  (8) 

Since  the  walker  spends  more  time  near  states  where  P  x  is 
large,  the  local  muximums  of  P  x  are  of  special  importance; 
they  are  called  the  stable  states  of  the  system. 

For  explicit  functions  a  ,  and  a  ,  we  can  easily  calculate 
the  function  P ,  on  a  digital  computer  directly  from  the  recur¬ 
sion  relation  equation  (6).  Figure  1  shows  a  computer¬ 
generated  plot  of  P  x  for  a  particular  model  in  which  equation 
(8)  is  a  cubic  with  three  distinct  roots.  This  system  evidently 
has  two  stable  states,  one  at  v ,  =  60  and  the  other  at.r-.  ~600. 
Although  it  is  not  evident  on  this  plot.  P  J.r)  is  strictly  positive 
for  all  .cat)  and  has  a  relative  minimum  at  m  =250. 


The  Fokker-Planck  Equation 

The  difficulties  encountered  in  trying  to  solve  the  master 
equation  have  prompted  some  investigators  to  seek  ap¬ 
proximating  equations  in  which  the  state  variable.*  is  regarded 
as  continuous  instead  of  discrete.  To  derive  such  an  ap¬ 
proximating  equation  for  equation  (3),  let  the  quantity 

'/;  |</  .(.*  +  I  )/»(.*  +  I  ,t)  +  a  (.*  -  1  )P(*  -1.0) 

be  added  and  subtracted  to  the  right  side.  A  simple  rearrange¬ 
ment  of  terms  then  yields 

HPU.nlH l  =  —(*•(.*+  1)P(.*+  I  ,f  )-*•(.*-  1  )P{X-  1  ,t)]l2 
+  '/i|*/(.*  +  I )  P(.v+  I . i )  -2 «(.*)  P(x,t) 

+atx—  I)  P(x-  !,/)],  (9) 

where 

*(-*)  =  a  Ax)-a  (.*)  and  of*)  =  aAx)+a  (x).  (10) 

Now,  if  x  were  a  continuous  variable,  the  two  terms  on  the 
right  side  of  equation  (9)  would  be  mutually  compatible  ap¬ 
proximations  of.  respectively,  the  first  jr-derivative  of 
-v(.*)P(.*,r)  and  the  second  .x-derivative  of  a(x)P(x.r): 
moreover,  these  approximations  should  be  good  provided  that 
r*>l.  We  thus  conclude  that  an  approximate  "continuized" 
version  of  the  master  equation  (3)  is 

£p(.*.r)  =  -  Jj(r(x)/>t*./)]  +  ^  |^j[at*)P(x.r)].  (**>!) 

(ID 

Equation  ( 1 1 )  is  called  the  Fokker-Planck  equation .  If  we 
multiply  through  by  x  and  then  integrate  over  x  by  parts,  we 
find  that  d<x  )Jdt=  <*■(*)),,  in  agreement  with  equation  (5). 
Mathematically,  equation  ( 1 1 )  is  identical  to  the  equation  gov¬ 
erning  a  one-dimensional  diffusion  process,  with  P(x.r)  play¬ 
ing  the  role  of  the  molecular  density  function,  vfr)  playing  (he¬ 
roic  of  the  average  molecular  drift  velocity,  and  a(x)!2  playing 
the  role  of  the  diffusion  coefficient.  In  the  particular  case  that 
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a  .(.v)  =  a  (x)=D,  equation  (II)  reduces  to  the  familiar  equa¬ 
tion  (iPj  dt  =D  if  PI  hx2 .  The  Fokker-Planck  equation  is  an 
exact  equation  for  certain  continuous-state  Markovian  pnx'es- 
ses;  however,  when  the  states  are  discrete  it  must  be  regarded 
as  an  approximation  that  is  good  only  for  single-step  processes 
with.v^l. 


The  Simulation  Algorithm 

In  many  ways,  the  function  P(x.t)  is  more  descriptive  of 
an  ensemble  of  walkers  than  of  an  individual  walker.  One  way 
of  focusing  on  an  individual  walker  is  to  simulate  the  random 
walk  on  a  computer.  The  key  to  doing  this  correctly  is  a 
function  // ,  defined  by 

p{T,fi\s,t)dT  =  the  probability  that,  given  the  walker  is  in 
state .v  at  time/,  it  will  step  next  to  state 
,v  +  /z  in  the  time  interval  (/  +  t,/  +  r  + 
dr).  (/^0,/x=  ±1)  (12) 

Technically  speaking,//  is  the  "joint  probability  density  func¬ 
tion"  of  the  random  variables  rand  //.  These  two  variables 
evidently  specify  when  and  to  where  the  walker's  next  step  will 
be. 

We  can  derive  a  formula  for//  in  the  following  way:  Since 
[</  ,(.v)r//+tf  UV//]  =</(.v)r//  is  the  probability  that  the  walker 
will  leave  its  current  slate  in  the  next  dr.  then  well  known 
probability  arguments  imply  that  exp|  -</(.v)r|  is  the  probabil¬ 
ity  that  the  walker,  in  state  \  at  time  /,  will  remain  there  during 
(/./  +  r).  Furthermore,  a  JxVlr  is  by  definition  the  probability 
that  the  w  alker,  in  state x  at  time  /  +  r,  w  ill  step  to  state x  +  fx  in 
the  next  dr.  The  probability  defined  in  equation  ( 12)  is  there¬ 
fore  just  the  product  exp|  trlxu^vk//.  Hence, 

//(  r,  /i|.v ,  / )  =  {//  „(.v  )/atx  )}{</(.»•  )exp( -</(.v )  r)  | 

I  />().//=  ±1).  (13) 


Equation  (13)  tells  us  (i)  that  r  and  //  are  mutually  inde¬ 
pendent  random  variables:  (ii)  that  r  is  exponentially  distrib¬ 
uted  with  decay  constant  </(.v);  and  (iii)  that  //  has  probability 
</.0)/</0)  of  being  equal  to  ±1.  It  is  relatively  simple  to 
"generate"  on  a  digital  computer  a  pair  of  random  numbers  t 
and  //according  to  these  specifications.  First,  we  generate  two 
random  numbers  r,  and  r.  from  the  unit-interval  uniform  dis¬ 
tribution;  then  we  take 


and 


t  =  |  l/«(.r)|log<  l/r,),  (14a) 


M  = 


|  -  I  if  a  U)>  r  Mix ), 
if  a  (.t)  s  r  ,</(.»  ). 


(14b) 


The  ability  to  properly  generate  a  random  pair  (r,//)  is 
evidently  all  that  we  need  to  simulate  the  random  walk.  Given 
initial  values  of  the  time  variable  /  and  the  state  variable  v,  we 
(i)  evaluate  the  functions  a  .(»)  and  a  (.< ).  (ii)  draw  two  ran¬ 
dom  numbers  r ,  and  r  .  from  the  unit-interval  uniform  distribu¬ 


tion,  (iii)  calculate  r  and  //  according  to  equations  ( 14a)  and 
( 14b),  and  (iv)  replace/  by  /  +  r and  x  by.v  +  /x.  Repeating  this 
cycle  of  steps  over  and  over  (on  a  computer),  and  periodically 
recording  the  instantaneous  state  of  the  walker,  we  obtain  a 
"realization"  of  the  random  walk  process,  if  the  uniform 
random  number  generator  is  "good,"  the  foregoing  simula¬ 
tion  algorithm  is  exact  with  respect  to  Equation  ( 1 );  moreover, 
it  is  easily  extendable  to  multivariable,  multistep  random 
walks. 

Figure  2  shows  the  results  of  eight  separate  simulations  of 
the  random  walk  whose  P .  function  is  plotted  in  Figure  I 
These  simulations  demonstrate  vividly  the  significance  of  the 
stable  states.  They  also  show  how  the  widths  of  the  stable  state 
peaks  in  thef*  curve  characterize  the  range  of  the  spontaneous 
fluctuations  of  the  walker  about  these  states. 

The  Functions  T+Xand  T  x 

We  let  t %  be  a  time  long  enough  after  the  random  walk 
starts  that  Ptx.l  J=P  Ax).  or 

PAx)=  the  probability  that  the  walker  will  be  in 

state  v  at  time  /,.  (15) 

Then  we  define  two  functions  T  \  and  T  „  by 

T  jtx)  *  the  average  time,  as  measured  from  /,. 

before  the  walker  first  executes  the  transition 
v— ».v±l.  (lb) 

We  can  obtain  an  explicit  formula  for  T ' ,(  v )  through  the 
following  line  of  reasoning:  equations  ( I )  and  (15)  imply  that 
/>,(.»)  x  a  .( x)dt  is  the  probability  that  the  walker  will  make  the 
transition  ,v— h  ±1  in  any  dt  interval  after  time/,.  This  implies 
that  the  time  interval  from  /x  to  the  first  occurrence  of  the 
transition  v—M  ±1  is  an  exponentially  distributed  random  vari¬ 
able  with  decay  constant  PAx)a  ,(.v).  Since  the  average  of  an 
exponentially  distributed  random  variable  is  just  the  reciprocal 
of  its  decay  constant,  we  conclude  that 

7'  ,(.v)  =  |«  .(.v)/\(.v)|  '.  (17) 

Equations  ( 17)  and  (b)  together  imply 

T  ,(.v)  =  r,(H),  (18) 

which  says  that  the  curve  of  7"  ,  is  just  the  curve  of  T shifted 
right  one  unit  of.v.  Il  follows  that,  for  most  practical  purposes, 
these  two  functions  can  be  regarded  as  one  and  the  same. 

If  we  multiply  equation  (b)  by  a.(x)  and  then  invoke 
equation  (17),  we  discover  that 

T  .(x)  =  |«  (.r)/«.(.r)J7'\(.t-  1).  (.v  =  1,2.  .  .  .)  (19) 

This  is  a  recursion  relation  for  T  .,  iteration  of  which  yields 

T  (ar)  =  7-  (0)11  [«  {n)/a .(«)].  0  =  1.2.  .  .  .)  (20) 

n  I 

The  quantity  T'j. 0)  can  be  calculated  simply  as  |c/.(0) 
PdO)]  !.  Comparing  equations  (19)  and  (20)  with  equations 
(b)  and  (7).  we  sec  that  the  functions  T\ are  just  as  easy  to 
calculate  as  the  function  Px. 


Then  /OS 
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The  "derivative"  of  T\  at.x  is,  using  equation  (19), 

T'M-T'J*- 1)  =  r  ^.r)t  1  -  a  Mia  (*>].  (21) 

We  see  from  this  that  the  extrema  of  T  \  are  those  points x  for 
which 

a  (x)  =  a  ,(.t  ),  (22) 

i.e. ,  those  points  for  which  the  walker  is  equally  likely  to  step 
right  or  left.  We  also  see  from  equation  (2 1 )  that  T • *  slopes  up 
or  down  at.t  according  to  whether  a  .(.v)  is  less  than  or  greater 
than  a  U )  That  implies  that  the  walker  always  has  a  bias  to 
move  "downhill"  on  the  7' \  curve;  moreover,  the  greater  the 
downward  slope,  the  stronger  the  bias.  The  functions  T  \  thus 
have  somewhat  the  character  of  a  potential  function  for  the 
random  walk.  In  particular,  it  follows  from  equations  (22). 
(8).  and  (17)  that  the  relative  minimums  of  T\  approximately 
coincide  with  the  relative  maximums  of  and  vice  versa. 
Therefore,  the  stable  states  of  the  system  can  just  as  well  be 
defined  as  the  relative  minimums  of  the  functions  7  • ,  and  T  ,. 

The  foregoing  points  are  illustrated  in  Figure  3,  which 
shows  a  plot  of  T  ,  for  the  same  system  considered  in  Figures  I 
and  2.  Because  of  equation  (18).  a  plot  of  T  „  would  kx>k 
virtually  identical.  The  T ' x  plot  shows  much  more  clearly  than 
the  P,  plot  the  true  nature  of  the  "barrier"  at  m  -250  that 
separates  the  stable  states  at  *,*=60  and  i :  =600 . 

First  Passages  and  Stable-State 
Transitions 

An  inspection  of  Figures  2  and  3  prompts  the  question  of 
how  long  it  would  take  a  walker  imtally  in  one  stable-state 
"valley"  to  beat  the  odds  and  make  it  over  the  "hill"  into  the 
other  stable-state  valley.  In  other  words,  what  is  the  average 
time  for  a  spontaneous  transition  from  one  stable  state  to 
another?  Such  a  question  could  arise  in  connection  with.  say. 
the  stability  of  an  electronic  memory  element  against  a  spuri¬ 
ous  change  of  stale,  or  the  stability  of  a  living  cell  against 
mutation. 

Let  us  pose  the  problem  in  a  slightly  more  general  way:  If 
the  walker  starts  in  any  state  x=u.  what  is  the  average  or 
expected  time  Hu—*v)  until  it  first  reaches  any  other  state  r  =v? 
We  need  not  require  u  and  v  to  be  stable  states,  but  for  the 
present  discussion  we  shall  require  that  v>u. 

To  calculate  r(t/— * »•>  we  first  introduce  some  auxiliary 
quantities,  all  defined  relative  to  the  walker's  “first  passage" 
from  u  to  »•:  t(x;u.v)  is  the  average  total  time  spent  in  state  v. 
n(x;u,v)  is  the  average  total  number  of  visits  to  state  .v,  and 
n .  (x.u.v )  is  the  average  total  number  of  x— ».v .  1  transitions.  In 
our  derivation  of  the  simulation  algorithm,  we  found  that 
l/«(.r)  is  the  average  time  the  walker  spends  in  state .r  on  any 
one  visit;  thus,  t(x;u.v)ln(x:u,v)=  l/«(.v).  We  also  found  that 
a .  ix)la(x)  is  the  probability  that  the  walker,  upon  leaving  state 
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Figure  3 

The  function  T\(x)  tor  the  system  of  figures  1  and  2  0 
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a  ,  will  jump  tostate.r±  1 ;  thus,  n . (a,«.v  )/rHx;u,v)=u  ,(x)la{x). 
Dividing  this  last  equation  into  the  preceding  one,  we  get 

Hx:u, v)  =  n  Ax:u. \  )lu  .(a).  (23) 

By  definition,  a  first  passage  from  ti  to  \  >u  must  have 
exactly  one  transition  i  -  I— *i ,  so 

n  .<!■-  I.n.r)  =  I .  (24a) 

Also,  the  number  of  a  -  I  — »a  transitions  minus  the  number  of 
a  -  I  *—x  transitions  must  be  exactly  zero  if  I  sa  su  and  exactly 
one  if  u  +  I  Sr  -  1 ;  thus, 

n  .(a  l.u.i)  -  n  (a v )  =  W(A-ttl. 

(l&vSf-l)  ( 24  b ) 

where  (Hz)  is  zero  if  ; s()  and  one  if  r'>0.  Dividing  equation 
(24a)  by  a  .0  -  I )  and  equation  (24b)  by  a  .(a  -  I ).  and  then 
making  use  of  equation  (23).  we  obtain 

/(»•-  l.d.t)  =  1  In  .0  I ).  (25a) 


/(»  l.n.i) 


a  (a  ) 
ti  .(a  -  I  ) 


/(A.d.r) 


«(A  -II) 

II  .(A  -  1  ) 


1  1  S|  1  ) 


(25b) 


Equations  (25a)  and  (25b)  constitute  a  complete  recursion  rela¬ 
tion  for  the  average  time  r(  v.«.r)  the  walker  spends  in  state  i . 
Once  these  times  have  been  calculated  for  all  >.  we  can  then 
compute  /(«—»»•)  as  their  sum: 

/(«—►!•)  =  1  Hx:ii. i  ).  (2b) 

Equations  (25a).  (25b).  and  (2b)  lend  themselves  nicely 
to  numerical  calculations  offdt— *rl  on  a  digital  computer  for 
explicit  functioned  .  and  a  :  in  that  way  .  one  can  show  ti  the 
model  in  Figures  1-3  that  the  average  time  for  the  stable  state 
transition  a,— »v  ,  is  m  1—»x;)  =  2.(W  x  10"  time  units.  Analyti¬ 
cal  calculation  of/(«— »r)  is  also  possible  w  ith  equations  (25a). 
(25b).  and  (2b).  For  example,  for  the  simple  case  u  Ax) 
=a  (a  )=D  considered  earlier,  one  can  easily  show  from  equa¬ 
tions  (25a)  and  (25b)  that  ttr-\:0.r)=  xll );  substituting  this 
into  equation  (2b)  gives  /<() — ■ >r)  =  r(r • +•  1  )/2D .  A  general 
analytical  iteration  of  equations  (25a  id  (25b)  is  algebrai¬ 
cally  messy  ;  however,  if  one  invokes  ti.  ulas  foi  P ,  and 
r  ,  in  equations  (7)  and  (17).  one  can  obta  : on  for 

rfv.d.i  )  that  allows  equation  (2b)  to  be  written  ;  ( -  ipact 
form 

tUi-*v)  =  I  (Pjf.)  -  7"‘jt«)  ).(«<»■)  (27) 

Finally,  we  may  note  that  equation  (27)  simplifies  even 
further  if  ii  and  i  are  the  stable  states  of  a  bistable  system. 
Referring  to  Figures  I  and  3.  and  putting  «-v,  and  i  =a  -  in 
equation  (27).  we  see  that  the  a -summation  contributes  sig¬ 
nificantly  only  for  a  under  the  a  ,-peak  in  the  P ,  curve,  while 
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ihc /(-summation  contributes  significantly  only  for/i  under  the 
m-peak  in  the  /  ,  curve.  Therefore.  the  average  time  for  a 
spontaneous  transition  from  one  stable  state  to  the  other  can  be 
well  approximated  by  the  product  of  (the  area  under  the  tuiluil 
stable  state  peak  in  the  P ,  curve}  times  {the  area  under  the 
intervening  barrier  peak  in  the  curve} — a  surprisingly  sim¬ 
ple  result. 

Concluding  Remarks 

Scientists  trained  in  classical  "deterministic"  modes  of 
analysis  are  often  ill-equipped  to  deal  with  the  uncertainties 
inherent  in  nondeterministie  systems,  f  or  these  systems, 
neither  the  relevant  dynamic  quantities  nor  the  fundamental 
equations  that  they  obey  are  obvious,  and  new  perspectives 
and  approaches  are  needed.  Simple  random  walk  models  like 
the  one  considered  here  are  important,  not  only  because  they 
approximate  some  real-world  systems,  but  also  because  they 
provide  a  context  in  which  new  perspectives  and  approaches 
can  be  developed  and  tested. 


Author 


Daniel  T  Gillespie  received  his 
H  A  trom  Kiee  University  in 
1460  ami  his  Ph  D.  (Physics) 
from  Johns  Hopkins  I'niversity 
in  I  46X  In  1 47  I .  alter  spending 
three  years  as  a  postdoctoral 
research  associate  al  the 
I'niversity  ol  Maryland's 
Institute  for  Molecular  Physics, 
he  joined  the  Research 
Department  ol  the  Naval 
Weapons  Center.  C  hula  l  ake. 
California,  as  a  research 
physicist.  He  is  now  head  ol  that 
department's  Applied 
Mathematics  Research  Branch 
His  current  research  interests 
are  in  kinetic  theory,  applied 
stochastic  process  theory,  and 
Monte  Carlo  methods 


Bibliography 

1.  Gillespie.  D  .  "A  General  Method  for  Numerically 
Simulating  the  Stochastic  Time  P.volution  ol  Coupled  Chemical 
Reactions."  Journal  of  ( omputalional  Wivun .  22  ( I 4761  p,  403. 

2.  Gillespie.  D  .  "A  Pedestrian  Approach  to  Transitions  and 
fluctuations  in  Simple  Noncquilihrium  Chemical  Systems." 
I’hwt  in.  Us  A  (  1474)  p  64 

Gillespie.  D  .  "Approximating  the  Master  liquation  by 
l-okker  Planck  Type  liquations  lor  Single- Variable  Chemical  Sys¬ 
tems."  Journal  ol  Chemical  Physics.  72  t  14X0)  p.  5367. 

4.  Gillespie.  D..  "A  Stochastic  Analysis  ol  the  Homogene¬ 
ous  Nucleation  of  Vapor  Condensation."  Journal  of  Chemical 
Phssn  \ .  74  1  I  US  | )  p.  66 1 . 

5.  Gillespie.  D..  "On  the  Calculation  of  Mean  first  Passage 
Tunes  for  Simple  Random  Walks."  Journal  ol  Chemical  Phssn  v . 
74  (14X1)  p.  5245 

6.  Gillespie.  D..  "  A  function  Csclul  for  Analyzing  Single- 
Step  One- Dimensional  Markov  Processes,”  Journal  of  Chemhal 
Physics.  76  (14X2)  p.  3762. 

7.  Calldatier.  R  .  "Stability  in  the  Dissipative  Steady  Stale.  " 
Physics  Taj  ax  (November  147X1  p.  2  1 . 

X.  Matheson.  I  .  D.  Walls  and  C.  Gardiner.  "Stoehastie 
Models  of  Pirst-Ordcr  Nonequi librium  Phase  Transitions  in  Chemi¬ 
cal  Reactions."  Journal  ol  Statistical  Physics,  12  (1475)  p  21. 

9.  Seshadri.  V  .  B.  West  and  K  l.indenberg.  "Analytic 
Theory  of  fxtrema  III  Results  for  Master  liquation  with  Applica¬ 
tion  to  l  ll l molecular  Decomposition,"  Journal  of  Cltciniial 
Phvsics ,  72  (14X0)  p.  I  145. 

10.  Weiss.  G  "first  Passage  Time  Problems  lor  One- 
Dimcnsional  Random  W  alks."  Journal  of  Statistical  Physics.  24 
1 14X1 1  p  5X7. 


52 


Vi/v  ill  Hi  ifiiri  h  Hi 


r 


1 


Research  Notes 


Healing  of  Propellant  Damage 

NWC  studies  of  the  damage-healing  phenomena  of  solid 
rocket  propellants  indicate  that  at  least  partial  healing  of  pro¬ 
pellant  damage  occurs  within  a  short  time  after  the  damage  is 
introduced.  The  results  of  these  studies  am  applicable  to  the 
rocket  motor  service-life  prediction  program  under  develop¬ 
ment. 

Composite  rocket  propellants  are  highly  filled  polymers 
'  ith  up  to  88'T  filler  consisting  mainly  of  oxidizers.  Any  finite 
deformation  of  a  propellant  can  cause  damage  to  its  micro¬ 
structure  in  the  form  of  molecular  chain  scissions  and  interfa¬ 
cial  debonding  or  dewetting.  A  damage  energy  criterion  de¬ 
veloped  at  NWC  was  applied  to  examine  the  healing 
phenomenon  of  propellant  damage.  (For  more  details,  see 
Journal  of  Applied  Polymer  Science .  Vol.  23.  Id7d.) 

Samples  of  both  unplasticized  and  plasticized  propellants 
were  predamaged  to  various  strains  and  stored  at  three  differ¬ 
ent  temperatures.  Several  sample  withdrawals  were  made  and 
the  remaining  damage  was  determined  as  a  function  of  time  up 
to  223  days.  Some  typical  results  are  illustrated  in  Figure  I . 


Figure  1 

Fractional  damage  energies  versus  storage  tune  at  three 
temperatures  of  plasticized  propellant  Mi  566  Numbers  m 
parentheses  indicate  number  ot  specimens  tested 
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The  healing  rate  or  decrease  in  damage  energy  is  highest 
during  the  first  few  weeks  after  damage  and  increases  with  the 
storage  temperature.  The  healing  then  ceases  for  a  longer  time 
period  and  even  reverses  at  lower  storage  temperatures  before 
it  finally  eontinues.  The  unplaslicizcd  propellant  showed  faster 
healing  properties  than  did  the  plasticized  propellant 

The  results  of  the  study  indicate  that,  while  at  least  partial 
healing  of  damage  occurs,  the  decrease  in  damage  energy  w  ith 
time  and  temperature  cannot  be  described  by  a  simple  kinetic 
law  over  a  longer  time  period.  The  breaking  of  molecular 
chains  is  a  major  contributor  to  damage.  The  free  radicals 
formed  are  reactive  species  and  are  expected  to  have  short 


lifetimes  because  of  the  presence  of  oxygen  and  other  ingre¬ 
dients  such  as  plasticizers,  inhibitors,  orchain-transfer  agents. 
Therefore,  the  conclusion  can  be  reached  that  the  healing  in¬ 
volves  a  number  of  chemical  reactions  that  result  in  secondary 
molecular  bonds,  which  then  form  new  primary  bonds  that  can 
support  a  mechanical  load. 

The  above  considerations  suggest  that  healing  lakes  place 
within  a  short  time  after  the  damage  is  introduced  into  the 
propellant.  The  experimental  results  confirm  rapid  changes  in 
the  propellant's  damage  level  at  the  beginning  of  the  healing 
period.  The  initial  modulus  recovers  l()0'3  as  the  damage 
energy  approaches  zero.  However,  the  ultimate  properties  of 
the  "healed"  propellant  are  changed — the  sires-,  at  failure  is 
slightly  increased  and  the  strain  at  break  reduced,  compared 
with  the  baseline  properties.  Further  property  changes  in  the 
same  direction  occur  as  the  propellant  undergoes  further 
aging. 

t Albert  H  l.cpic.  Vtf'O 

Analytical  Chemistry  and  Pollution 
Abatement 

Over  the  years  Naval  Weapons  Center  (NWO  chemists 
have  synthesized,  analyzed,  and  studied  a  number  of  new 
high-energy  chemical  compounds  for  use  in  propellants,  ex¬ 
plosives.  and  pyrotechnics.  One  recent  research  effort  led  to 
the  production  of  a  small  analytical  kit  to  help  the  Department 
of  Defense  monitor  the  pollution  levels  ot  trinitrotoluene 
(TNT)  and  ammonium  picrate(AP)  at  naval  bases  where  these 
explosives  were  once  dumped.  These  explosive  compounds 
can  seep  into  water  tables  and  appear  in  well  w  ater  w  here  even 
It)  pans  per  billion  is  toxic. 

The  chemical  pollution-detection  method  developed  at 
NWC  depends  upon  complex  chemistry,  but  the  final  package 
is  a  small  kit  that  can  be  used  by  an  engineer,  technician,  or 
sailor  who  knows  nothing  about  chemistry.  The  kit  contains  a 
syringe  pump  and  glass  tubes  filled  w  ith  the  detector  chemi¬ 
cals.  A  given  volume  of  water  is  pumped  through  the  indicator 
tube,  and  a  color  appears  as  the  w  ater  passes  through  the  tube 
The  length  of  the  color  section  measures  the  concentration  of 
TNT  in  the  water.  The  tubes  arc  expendable,  and  tubes  of 
different  diameter  are  used  for  different  concentration  ranges. 
Doubling  the  given  amount  of  water  doubles  the  color  length 
for  a  given  concentration.  A  solution  of  know  n  TNT  concen¬ 
tration  is  included  in  the  kit  for  calibration  of  the  length  of  the 
colored  section  obtained  versus  the  concentration  of  the  know  n 
solution. 

The  chemistry  is  interesting  and  depends  upon  the  forma¬ 
tion  of  Meisenheimer  anions  These  colored  anions  arc  then 
trapped  at  cationic  sites  on  the  surface  of  ion-exchange  resins 
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The  structure  of  Meisenheimer  anions  is  not  completely  cer¬ 
tain  .  hut  present  evidence  w  ould  suggest  the  following  chemis¬ 
try: 


NO,  NO., 

2  2  (I) 


This  structure  is  now  favored  over  a  structure  w  ith  one  proton 
removed  by  the  hydroxide  anion.  Evidence  from  NWC  work 
suggests  that  the  anions  may  form  dimers  and  larger  poly  liters 
that  block  the  resin  pores  Various  si/e  polymers  may  account 
for  color  variations  from  red  to  purple. 

The  resin  is  chemically  a  polystyrene  chain  with  side 
groups  of  phenyl  rings  having  tertiary  amine  groups.  This 
chemistry  is  depicted  as: 


The  chloride  anion  is  replaced  by  the  more  firmly  bound 
Meisenheimer  anion  In  the  case  of  AP.  it  is  the  y  ellow  picrate 
anion  that  displaces  the  chloride  anion.  NWC  researchers  are 
study  mg  methods  to  alter  and  intensify  the  color  of  the  picrate 
ion 

The  chemical  reaction  1 1 )  is  carried  out  in  a  short  section 
ot  the  tube  filled  with  a  mixed  alkaline-earth  oxide  section. 
Here  the  pH  is  high  enough  to  drive  all  of  the  TNT  to  form 
Meisenheimer  anions  ( Irigmally .  this  section  contained  only 
Cat),  but  the  ionic  strength  ot  NaCl  in  brackish  water  pre¬ 
vented  the  formation  of  Meisenheimer  anions.  Now  ,  with  the 
addition  ot  Cat).  Sr(>.  and  Ba().  the  tubes  can  be  operated  in 
areas  such  as  the  tidal  Hals  of  Virginia. 

f-.nv  ironmental  engineers  have  long  used  color  indication 
tubes  for  detecting  air  pollutants,  but  tubes  for  detecting  water 
pollutants  are  a  new  tool  in  their  arsenal.  With  this  basic  idea 
now  established,  it  is  likely  that  similar  tubes  to  detect  other 
pollutants  will  be  developed  for  DOD. 
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Organometallic  Precursors  for  the 
Production  of  Mercury  and  Cadmium 
Tellurides 

Researchers  at  the  Naval  Weapons  Center  (NWC)  are 
investigating  the  use  of  organometallic  precursors  for  the  pro¬ 
duction  of  ll-VI  semiconductors  in  general  and  for  mercury 
and  cadmium  tellurides  specifically.  The  use  of  these  or- 
ganometallic  precursors  is  being  pursued  for  the  preparation  of 
thm-film  semiconductors  with  the  desired  stoichiometry  (ratio 
of  numbers  of  atoms  of  mercury,  cadmium,  and  tellurium  to 
each  other)  and  purity  over  large  areas. 

Infrared  sensors  constructed  of  thin  crystalline  films  of 
mixed  mercury -cadmium  telluride  deposited  on  a  substrate  of 
single-crystal  cadmium  telluride  have  the  advantage  of  being 
able  to  operate  in  two  different  atmospheric  w  indows.  The  thin 
films  must  be  of  high  chemical  purity  and  must  be  applied  to 
substrates  as  free  from  crystalline  defects  as  possible:  the 
stoichiometry  of  the  films  must  be  closely  controlled,  and  the 
film  must  he  uniform  over  the  entire  surface.  Further,  the  films 
may  need  to  be  deposited  over  larger  and  larger  areas  as  new 
requirements  for  infrared  sensors  and  other  semiconductor- 
film  devices  arise. 

To  meet  requirements  for  high  purity  and  closely  control¬ 
led  stoichiometry  over  large  areas.  NWC  is  investigating  and 
developing  techniques  using  organometallic  precursors  to  pre¬ 
pare  semiconductor  compounds. 

The  first  step  in  the  preparation  of  semiconductors  using 
organometallic  precursors  involves  the  synthesis  and  charac¬ 
terization  of  specific  model  compounds.  The  model  com¬ 
pounds  of  immediate  interest  are  represented  as  RMXR’  (M  = 
/. n.  Cd.  or  Hg;  \  =  S.  Se.  orTe;  R  and  R'  =  organic  groups). 
The  use  of  these  compounds  is  a  unique  approach  to  the  prep¬ 
aration  of  ll-VI  materials;  they  have  the  required  stoichiometry 
and  a  metal-chalcogenide  bond  has  already  been  formed.  In 
fact,  they  arc  similar  to  intermediates  postulated  during 
chemical-vapor  deposition  of  semiconductors. Compounds  of 
this  type  may  be  purified  chemically  to  almost  any  desired 
level  of  purity  . 

NWC  chemists  have  prepared  CaHsZnSCiHu. 
C/HvCdSCiH*.  <C«H.-,SehHg.  C.iH.vHgSeC.iH-,.  and 
(C«H.vTe)jHg.  and  have  also  attempted  to  synthesize 
C«H.-,HgTeC«H.-,  and  C«.H.-,CdTeC«Hs.  Although' the  latter 
two  compounds  have  not  been  fully  characterized .  compelling 
evidence  exists  that  the  attempts  have  been  successful. 

Another  aspect  of  the  Center's  work  is  the  use  of  light 
(photolysis)  or  of  temperatures  below  100°C  to  cause  decom¬ 
position  of  the  model  compounds  and  formation  of  the  ll-VI 
semiconductor  material.  The  use  of  these  relatively  mild  con¬ 
ditions  (temperatures  as  high  as  900°C  are  employed  in  other 
decomposition  methods)  greatly  reduces  the  migration  of  ions 
in  the  crystalline  materials  and  should  he  a  useful  technique 
when  devices  of  these  materials  are  manufactured.  The  forma¬ 
tion  of  HgSc.  HgTc.  and  CdTc  under  the  mild  conditions 
described  has  been  observed  at  NWC. 
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The  HgSe  was  prepared  by  photolysis  of 
CsHsHgSeCsHs,  dissolved  in  a  nonaqueous  solvent,  at  room 
temperature.  A  broadband  ultraviolet  source  was  used,  and 
oxygen  was  rigorously  excluded.  Under  the  same  conditions 
(C«HsSehHg  and  (CnHsTeViHg  showed  no  photolytic  reac¬ 
tion,  although  in  the  presence  of  sunlight  and  oxygen  from  the 
air  the  (CuHr.TehHg  reacted  to  give  (C«Hs)aHg  and  some 
complex  organic  tellurium  oxides. 

The  HgTe  and  CdTe  were  prepared  from  the  putative 
precursors  C«H.-,HgTcC«H.-,  and  C«H.-,CdTeC«H.-,.  The  HgTe 
was  prepared  by  heating  benzene  or  tetrahydrofuran  solutions 
to  about  80°C,  and  also  by  photolysis  at  about  30°C  under  the 
previously  mentioned  conditions.  The  CdTe  was  prepared  by 
decomposition  at  room  temperature. 

(Wayne  E.  Thun.  NWC) 

The  Role  of  Energetic  Polymers  in 
Propellants  and  Explosives 

There  is  an  ever-increasing  need  to  enhance  the  perfor¬ 
mance  of  energetic  materials,  particularly  propellants  and  ex¬ 
plosives.  in  both  military  and  non-military  applications.  Typi¬ 
cal  rocket  propellants  and  plastic-bonded  explosives  (PBX) 
contain  binders  of  inert  polymers.  While  the  inert  binders 
desensitize  the  hazardous  ingredients  and  provide  fluidity  (dur¬ 
ing  processing),  elongation,  and  tensile  strength,  they  also 
degrade  energy  release.  Although  efforts  are  being  made  to 
develop  superior  ingredients,  stable  and  compatible  com¬ 
pounds  that  have  appreciably  more  energy  than  the  current 
ingredients  of  ammonium  perchlorate,  eyclotetrumethylene 
tetramtraimne  (HMX),  and  aluminum  (Al)  are  unlikely  to  be 
developed  in  the  near  future.  Therefore,  the  use  of  energetic- 
binders  may  be  the  most  logical  approach  to  the  attainment  ol 
higher  performance. 

A  family  of  energetic  binders  (azido  polymers)  with 
polyether  backbone  has  been  developed  that  increases  the  per¬ 
formance  of  propellants  and  explosives  without  degrading 
safety ,  processing,  and  mechanical  properties.  The  binders  arc 
easily  prepared  in  the  form  of  pourable  liquids  in  the  2000  to 
4000  molecular  weight  range.  Polyethers  have  desirable 
rheological  properties  and  tend  to  be  soft  and  rubbery,  a  con¬ 
sequence  of  the  tendency  for  free  rotation  of  the  ether  oxygen 
to  carbon  linkage. 

Energetic  polymers  offer  the  potential  for  increasing  the 
burning  rate  of  normally  low-buming  propellants  containing 
HMX.  Mechanical  and  safety  properties  can  be  enhanced  be¬ 
cause  higher  volume  fractions  of  polymers  can  be  formulated 
without  degrading  performance,  a  consequence  of  the  in¬ 
creased  energy  content  of  the  polymer.  Studies  of  rocket  prop¬ 
ellants  of  equal  performance  indicate  that  the  safety  properties 
are  improved  as  the  level  of  HMX  is  reduced  and  the  energy 
contribution  of  the  binder  is  increased  by  the  addition  of  a  more 
energetic  plasticizer. 

The  results  of  recent  research  at  NWC  arc  summarized  in 
the  following  paragraphs. 


Performance 

When  inert  binders  are  replaced  by  denser  energetic 
polymers  in  propellant  compositions,  performance  is  en¬ 
hanced  because  of  the  gain  in  enthalpy  release.  Since  density 
as  well  as  enthalpy  release  are  involved  in  the  calculated  deto¬ 
nation  pressure,  the  denser  energetic  poly  mers  produce  a  rela¬ 
tively  greater  performance  increase  in  explosive  composition 
than  they  do  in  rocket  propellants 

Ballistic  Properties 

Energetic  polymers  increase  the  burning  rate  of  propel¬ 
lants  and  reduce  the  extent  of  aluminum  agglomeration.  Burn¬ 
ing  rates  are  generally  enhanced  as  the  azido  content  of  the 
polymer  increases 

Chemical  and  Thermal  Stability 

The  azido  polymers  seem  to  be  compatible  with  the  fol¬ 
lowing  typical  propellant  ingredients  as:  energetic  plasticizers 
containing  mtro.  fluorodinitm.  and  nitrate  ester  groups; 
oxidizers  such  as  ammonium  nitrate  and  ammonium  perchlo¬ 
rate  (AP);  mtramines  such  as  RDX  (a  powerful  but  sensitive- 
explosive)  and  HMX;  burning  rate  depressants  such  as 
oxamide.  and  fuels  such  as  the  thermogenic  light  metals  Al. 
Mg.  etc.  The  thermal  stability  of  the  prepolymers,  cured 
polymers,  and  propellant  and  explosive  compositions  was 
found  to  be  superior  to  that  of  nitrate  esters. 

Mechanical  Properties 

The  mechanical  properties  of  a  variety  of  compositions 
containing  energetic  polymers,  plasticizers.  HMX.  AP.  and 
Al  were  determined;  and  these  compositions  were  found  to 
have  excellent  elongation  and  adequate  tensile  properties 
Low  glass  transition  temperatures  were  exhibited,  an  indica¬ 
tion  that  the  propellants  will  show  gixxl  mechanical  properties 
at  low  temperatures. 

Conclusions 

The  current  results  of  the  continuing  investigation  have 
verified  the  predictions  that  the  use  of  energetic  binders  in 
propellants  and  explosives  will  ( 1 )  increase  energy  release,  ( 2) 
permit  ballistic  tailorability ,  (3)  enhance  safety  properties,  and 
(4)  improve  the  mechanical  properties  of  PBXs.  There  is 
growth  potential  for  propellants  and  explosives,  as  well  as  the 
potential  for  improved  energy  release  in  the  areas  of  gun  pro¬ 
pellants  and  pyrotechnics. 

( Russell  Reed  Jr..  NWC.  Max  L.  Chan,  NWC.  and  Richard  S. 
Miller.  ONR) 


Three-  INS* 


55 


Research  Notes 


Radar  Target  Imaging  and  Inverse  Vector 
Scattering 

The  signal  processing  techniques  for  the  detection  and 
tracking  of  radar  aircraft  echoes  have  been  thoroughly  investi¬ 
gated  over  the  years,  and  the  optimal  methods  are  well  known. 
However,  the  signal  processing  techniques  for  the  classifica¬ 
tion  (friend,  foe,  or  neutral)  and  identification  (for  example. 
F-4  or  MIG-23)  of  aircraft  targets  are  still  in  the  early  stages  of 
research.  The  Naval  Weapons  Center  (NWC)  is  investigating  a 
classification  and  identification  technique  that  uses  the  polari¬ 
zation  characteristics  of  the  radar  echo,  along  with  a  modest 
degree  of  Doppler  processing,  to  form  an  image  of  the  convex 
portions  of  the  aircraft  (such  as  the  fuselage  or  engine  pods). 
Then,  by  determining  the  diameter  and  length  of  the  image  (or 
perhaps  by  extracting  other  measurements  from  the  image), 
the  type  of  aircraft  can  be  identified. 

This  technique  employs  a  direct  attack  on  the  solution  of 
the  inverse  scattering  problem  by  using  the  vector  nature 
(polarization)  of  the  echo  signal.  This  method  is  in  contrast  to 
most  current  target  identification  methods,  which  require  a 
vast  file  of  experimentally  taken  radar  signals  on  actual  targets. 
That  file  is  then  searched  for  the  best  match  with  an  unknown 
radar  echo.  No  such  file  is  required  for  the  NWC  technique. 


The  depolarization  of  a  radar  signal  by  a  scattering  sur¬ 
face  provides  a  means  by  which  the  local  curvature  properties 
of  the  surface  may  be  determined.  Furthermore,  a  classic  prob¬ 
lem  in  differential  geometry  (Christoffel-Hurwitz)  deals  with 
the  reconstruction  of  such  a  surface  from  a  knowledge  of  this 
kind  of  information,  and  a  differential  equation  relating  these 
local  measurements  to  the  surface  has  long  been  established. 
Unfortunately,  its  exact  solution  requires  an  exhaustive 
amount  of  input  data,  more  than  could  be  realized  by  any 
expected  radar  encounter. 

It  is  possible,  however,  to  attack  not  the  differential  equa¬ 
tion  but  the  corresponding  variational  problem  and  thereby 
allow  for  a  solution  on  a  finite  data  set.  In  this  way.  a  solution 
that  best  fits  the  known  data  can  be  found. 

A  FORTRAN  code  employing  this  “finite-element"  in¬ 
version  algorithm  has  been  constructed  and  tested  on  synthetic 
data.  The  effect  of  perturbed  or  corrupted  data  on  this  al¬ 
gorithm  has  also  been  investigated,  as  well  as  the  extent  to 
which  the  data  set  can  be  limited. 

Future  plans  call  for  evaluating  the  technique  first  on  a 
94-GHz  experimental  radar  at  NWC  using  scale  model  targets, 
and  then  using  a  suitable  polari metric  radar  on  full-size  aircraft 
targets. 

( Brett  Borden  and  Robert  J.  Dinger,  NWC) 
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The  Naval  Weapons  Center  (NWC)  at  China  Lake.  lo¬ 
cated  1 50  miles  north  of  Los  Angeles,  covers  a  million  acres  of 
California's  Mojave  Desert.  In  1943,  a  Navy  and  California 
Institute  of  Technology  team  created  a  thriving  scientific  facil¬ 
ity,  known  as  the  Naval  Ordnance  Test  Station  (NOTS).  in  this 
vast  and  sparsely  populated  region.  Here  the  Navy  success¬ 
fully  met  the  urgent  wartime  need  for  aircraft  ordnance  de¬ 
velopment  and  testing  The  early  partnership  established  a 
pattern  of  cooperative  interaction  between  civilian  scientists 
and  experienced  military  personnel  that,  in  the  ensuing  four 
decades,  has  made  NWC  one  of  the  preeminent  research  in¬ 
stitutions  in  the  country. 

Toward  the  end  of  the  1940s,  NOTS  began  research  on  a 
fire-control  system  that  in  time  became  Sidewinder,  the  na¬ 
tion's  first  passive  infrared-homing  air-to-air  guided  missile. 
With  the  advent  of  the  Korean  conflict  in  1950,  NOTS  rapidly 
gained  cognizance  over  an  extensive  catalog  of  air-to-air  rock¬ 
ets,  guided  missiles,  air-launched  torpedoes,  and  an  array  of 
guns,  bombs,  and  fuzes.  By  the  late  1950s.  research  expanded 
into  such  diverse  fields  as  antisubmarine  weaponry  and 
weather  modification. 


American  involvement  in  Vietnam  in  the  1960s  quick¬ 
ened  the  tempo  of  activities  at  NOTS.  and  a  new  generation  of 
"smart"  bombs,  cluster  weapons,  and  night  attack  systems 
was  developed.  Electronic  warfare  received  major  attention, 
and  the  Station  also  made  important  technical  contributions  to 
Navy  electronic  countermeasures  and  counter-countermea¬ 
sures  capabilities. 

In  1967.  NOTS  and  the  Naval  Ordnance  Laboratory. 
Corona,  joined  to  form  the  Naval  Weapons  Center.  The  Center 
acquired  the  National  Parachute  Test  Range  (NPTR)  function 
in  1979  and  now  serves  as  a  national  parachute  test-and- 
evaluation  facility. 

Today,  NWC  carries  out  the  complete  development 
process  -from  basic  research  to  prototype  hardware  to  test  and 
documentation — on  over  800  separate  programs  or  tasks.  The 
Center's  major  programs  include:  Sidewinder,  Sparrow,  and 
Phoenix  air-to-air  missiles:  Harpoon  and  Walleye  air-to- 
surface  weapons;  the  high-speed  antiradiation  missile 
(HARM);  fuel-tank  recovery  parachutes  for  the  space  shuttle; 
and  avionics  support  for  a  variety  of  tactical  aircraft. 


